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This paper presents some new methods of computing the sedimentation radius. The methods 
are founded on a curve for spheres, experimentally determined and constructed on the basis 
of the coefficient of resistance as a function of Reynolds number. When the validity of Stokes’ 
law ceases, the sedimentation’ radius obtained by Stokes’ formula is graphically corrected. 
The amount of deviation of Stokes’ law from experimental data of higher Reynolds numbers 
has been statistically determined. Some new formulas for computing the resistance and the 
radius of a sphere of a given terminal settling velocity are presented. 





SoME NEW TERMS AND NOTATIONS 


_proead introduced in 1868 the term hy- 
draulic value (hydraulischer Werth), meaning 
thereby the diameter size of a quartz sphere 
having the same settling velocity as a given 
particle in water as sedimentation medivm. Thus 
the hydraulic value has no bearing on the actual 
size of the particle. An improvement was made 
by Odén? in 1915 by introducing the conception 
equivalent radius, i.e., the radius of a sphere of 
the same material and of the same settling 
velocity as a given particle. The sedimentation 
medium is not expressively stated, but one 
infers from Odén’s writing that water is used as 
sedimentation fluid. Although I believe that 
Odén by “‘the same material’’ means the same 
Specific gravity, the definition is not sharp. The 
specific gravity value, but not the material, 
enters the sedimentation formula. The term 
equivalent radius as synonymous with the radius 
of a sphere of equivalent settling velocity is 


1 Em. Schéne, Ueber einen neuen A pparat fiir die Schlim- 
manaiyse, Zeits. f. analyt. Chemie 7, 29 (1868). 

2S. Odén, Eine neue Methode zur mechanischen Boden- 
analyse, Int. Mitt. f. Bodenkunde 5, 257 (1915). 


slightly confusing, if one wishes to use the 
expressions equivalent volume, equivalent surface, 
etc. In conclusion, I believe that a new term 
sedimentation radius is justified. Definition: The 
sedimentation radius is the radius of a sphere of 
the same specific gravity and of the same terminal 
uniform settling velocity as a given particle in the 
same sedimentation fluid. The term does not 
imply any special kind of sedimentation medium, 
and the kind of fluid used must be stated in 
each case. The values obtained for the sedimenta- 
tion radius in different fluids are not com- 
parable above Reynolds number 0.2, except when 
the particles are perfect spheres. 

The diameter and the radius of a sphere of 
the same volume as a given, spherical or non- 
spherical, solid have been previously introduced 
as the true nominal diameter and the true nominal 
radius, respectively.*: 4 § 


3H. Wadell, Volume, Shape and Roundness of Rock 
Particles, J. Geology 40, 443 (1932). 

4H. Wadell, Sphericity and Roundness of Rock Particles, 
J. Geology 41, 310 (1933). 

5H. Wadell, The Coefficient of Resistance as a Function 
of Reynolds Number for Solids of Various Shapes, J. Frank. 
Inst. 217, 459 (1934). 
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The actual radius is an emphasized expression 
for the radius of a perfect sphere. The term is 
used as a contrast to the Stokes radius, which is 
the radius of a perfect sphere computed according 
to Stokes’ law. The Stokes radius is always 
smaller than the actual radius for Reynolds 
numbers higher than 0.2. The Stokes sedimenta- 
tion radius signifies that the sedimentation radius 
has been computed according to Stokes’ law. 
The Stokes velocity signifies the terminal uniform 
settling velocity computed according to Stokes’ 
formula for spheres. 

The conception degree of true sphericity has 
been previously introduced’: *:° as a sedimento- 
logical shape expression based on the isoperi- 
metric property of a sphere and expressed by 
the formula: ~=s/S, where y is the degree of 
true sphericity, s the surface area of a sphere of 
the same volume as the solid, and S the actual 
surface area of the solid. The maximum value 
obtained by this formula is 1.000, which is the 
degree of true sphericity of a sphere. 

The easiest shape position is that position in 
which a solid of a given shape attains its 
maximum terminal uniform settling velocity. 
The fardiest shape position signifies the position 
of minimum velocity. The average shape position 
is understood to be a great number of positions 
in which the solid attains different velocities, 
the mean of which constitutes the average velocity. 
Except for spheres, all shapes have an easiest and 
a tardiest shape position. 


THE COEFFICIENT OF RESISTANCE AS A FUNCTION 
OF REYNOLDS NUMBER FOR SOLIDS 
OF VARIOUS SHAPES 


The force, acting upon a falling particle 
immersed in a fluid, equals the effective weight 
of the solid. When the particle has attained its 
terminal velocity this force (F) is exactly bal- 
anced by the fluid resistance, R: 


F=R=Mg=(4/3)cr,3(ps— pr)g, (1) 


where is the mass, 7, the true nominal radius, 
ps the specific gravity of the solid, pr the 
specific gravity of the fluid, and g the acceleration 
of gravity. It follows that the resistance, R, 
must always be the same for solids of the same 
specific gravity, of the same true nominal radius 
(or the same volume) and immersed in the same 
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fluid, irrespective of shape and shape position 
of the solids, nota bene, when the terminal 
velocities of respective solids have been at- 
tained. The terminal velocity, which also is the 
maximum velocity, varies with the shape and 
shape position and also with the Reynolds 
number. This variation is conveniently expressed 
by the so-called coefficient of resistance. A slight 
change, by inserting the true nominal radius in 
conventional formulas for the coefficient of re- 
sistance and the Reynolds number for spheres, 
has been previously suggested. The formulas 
become thereby applicable to all shapes and 
shape positions, and the coefficient of resistance 
as a function of Reynolds number expresses the 
influence which the properties of the fluid and 
of the solid have on the terminal uniform 
settling velocity. The formulas 1ead: 


Cr=8r,.(ps— pr)g/3v’ pr, (2) 
R=2r,vpr/u=2r,v/v, (3) 


where Cpr is the coefficient of resistance, 7, the 
true nominal radius, ps the specific gravity of 
the solid, pr the specific gravity of the fluid, u the 
viscosity coefficient of the fluid, v the kinematic 
viscosity u/pr, g the acceleration of gravity, and 
R Reynolds number. 

The coefficient of resistance as a function of 
Reynolds number for spheres is shown in Figs. 
1, 2,3 and 4. The curve represents the results ob- 
tained by Allen,® Arnold,’ Liebster,* Schmiedel,’ 
Wieselsberger, Betz and Prandtl.!° Ladenburg’s" 
correction for the influence of the wall of the 
sedimentation vessel has been taken into con- 
sideration, whenever the corrected experimental 
data have been available. Schiller’s'® tables have 


6H.S. Allen, The Motion of a Sphere in a Viscous Fluid, 
Phil. Mag. [5] 50, 323, 519 (1900). 

7H. D. Arnold, Limitations Imposed by Slip and Inertia 
Terms upon Stokes’ Law for the Motion of Spheres Through 
Liquids, Phil. Mag. [6] 22, 755 (1911). 

8H. Liebster, Uber den Widerstand von Kugeln, Ann. d. 
Physik [4] 82, 541 (1927). 

9 J. Schmiedel, Experimentelle Untersuchungen iiber die 
Fallbewegung von Kugeln und Scheiben in reibenden Fliis- 
sigkeiten, Phys. Zeits. 29, 593 (1928). 

10°C, Wieselsberger, A. Betz, L. Prandtl, Ergebnisse der 
——" Versuchsanstalt zu Géttingen, Part II 

1923). 

1 R. Ladenburg, Uber den Einfluss von Wéanden auf die 
Bewegung einer Kugel in einer reibenden Flissigkeit, Ann. d. 
Physik [4] 23, 447 (1907). 

LL. Schiller, Fallversuche mit Kugeln und Scheiben, 
Handb. d. exp. Physik (Ed. by W. Wien and F. H. Harms), 
4, No. 2, Hydro- und Aerodynamik, pp. 337-387 (1932). 
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Fic. 1. Coefficient of resistance as a function of Reynolds number for spheres. 
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Fic. 2. Coefficient of resistance as a function of Reynolds number for spheres. 


also been consulted. The dotted curves represent 
the coefficient of resistance as a function of 
Reynolds number [Cr=f(R)] according to 
Stokes’ and Oseen’s formulas for spheres. Chart 
Fig. 5 illustrates the dispersion of the experi- 
mental data, expressed by the Cr =f(R) accord- 
ing to results obtained by Allen,* Liebster,*® 
Schmiedel,® Wieselsberger'!® and Lunnon." All 
data are represented without Ladenburg’s or 


1R. G. Lunnon, Fluid Resistance to Moving Spheres, 
Proc. Roy. Soc. (London) A118, 680 (1928). 


Faxén’s" corrections. Allen’s results refer to his 
spheres of amber and of steel in water. The 
data for the single dot representing Lunnon’s 
result have been taken from a table by Schiller." 

It has been shown in a previous paper® that 
the curves, constructed on the basis of the 
Cr=f(R) for solids of various shapes (sphericity 
values), become asymptotic towards lower Rey- 
nolds numbers. Within the region of laminar 

44H. Faxén, Die Bewegung einer starren Kugel lings der 


Achse eines mit 2dher Fliissigkeit gefiillten Rohres, Arkiv f. 
Matematik, Astronomi och Fysik 17, No. 27 (1923). 
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motion of the fluid about the solid, i.e., within 
the region of Stokes’ law, the Ce=/f(R) for a 
given shape and shape position appears to take 
a linear trend parallel with the linear function of 
Stokes’ law for spheres. On account of rotational 
motion of the fluid the effect of the shape in- 
creases with increasing Reynolds number, so 
that the distance between the curve for spheres 
and the Cre=/(R) for a given shape and shape 
position increases with increasing R-value. As a 
consequence, the difference in size between the 
sedimentation radius, 7,, and that of the true 
nominal radius, 7,, increases with increasing 
Reynolds number for a given shape and shape 
position. An example is given in Table I. Solids 
Nos. 2, 3 and 4 have the same shape positions 
and almost the same shape (y¥ = 0.2182 to 0.2375). 
The Reynolds number increases from 0.1389 
(solid No. 2) to 7881 (solid No. 4). 1? Among 
the symbols in the table th represents the thick- 
ness of a disk, rg the radius of the flat face of the 
disk, r, the true nominal radius, 7,5; the Stokes’ 


%V. Pernolet, A l'étude des préparations mécaniques des 
minerais, ou expériences propres a établir la théorie des 
différent systémes usités ou possibles, Ann. des Mines [4] 20, 
390 (1851). 

16H. Wadell, Shape Determinations of Large Sedimental 
Rock- Fragments, Pan-American Geologist 61, 187 (1934). 
Unfortunately several misprints occur in this paper. The 
integration symbol, /, is used throughout the paper 
instead of f as symbol for function. Eq. (8) on page 203 
should read: ¥ =d,/Ds5+0.1. 
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Fic. 3. Coefficient of resistance as a function of Reynolds number for spheres. 
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sedimentation radius, and 7, the sedimentation 
radius determined by methods described in a 
latter part of this paper. For the solids Nos. 2, 
3 and 4 the ratios 7,/r, and volume of a sphere of 
radius r,/volume of a sphere of radius r, are as 
follows: 





Solid Tn Volume of sphere of radius r,, 
No. rs Volume of sphere of radius r, 
2 1.6 4.6 
3 2.4 16.0 
4 12.3 1870.5 
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Fic. 4. Coefficient of resistance as a function of Reynolds 
number for spheres. 
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100 
Allen 
Wieselsberger 
Liebster 
Schmiedel/ 
10 Lunnon 
10 
01 
Of 1.0 10 100 1000 10000 
Fic. 5. 
TABLE I. 
"sSt 
Geometric or Volume v 
Solid shape th. Yd Yn rs cc Mate- cm/ 
No. classification v cm cm cm cm ((4/3) ry’) rial Ps PP “ v sec. R CR 
1 Circular disk 0.1219 0.02068 0.99900 0.24922 0.07176 0.064837 Steel 7.674 1.199 0.4417 0.3684 11.55 15.627 26.393 
(Schmiedel 13) 
2 Circular — 0.2182 0.01050 0.20255 0.06862 0.04127 0.0013533 Steel 7.674 1.223 1.7012 1.391 1.408 0.1389 477.60 
(Schmiedel 1) 
3 Circular disk 0.2182 0.010497 0.20250 0.06860 ©.02843 0.0013523 Steel 7.674 1.150 — 0.0906 7.569 11.46 17.77 
(Schmiedel 10) 
4 Circular disk 0.2375 0.125 2.100 0.74495 0.060463 1.7318 Lead 11.150 0.9986 — 0.0106 56.2 7881. 6.269 


(Pernolet 5) 








All values above refer to the tardiest shape 
position, i.e., with the flat face of the disk 
normal to the direction of translation. Table I 
shows that rg>rn>r7r;s. It is interesting to note 
that the difference between the true nominal 
radius and the sedimentation radius is com- 
paratively slight at low Reynolds numbers. Thus 
the actual size of a very flat particle, even in 
the tardiest position, can be determined with a 
fair accuracy by the sedimentation method, 
when Stokes’ law is applicable, i.e., below Rey- 
nolds number 0.2. If the effective weight of a 
particle in water as sedimentation medium is so 
high that Stokes’ formula is not applicable, the 
sedimentation can be brought within the validity 
of Stokes’ law, either by employing heavy liquids, 
thereby reducing the effective weight, the settling 


velocity and the Reynolds number, or else by 
using liquids of high viscosity with practically 
the same result, namely, reduced velocity and 
lower Reynolds number. These methods apply 
to investigations, the aim of which is to obtain 
a size determination of highest possible accuracy 
for non-spherical particles. 


COMPUTATION OF THE SEDIMENTATION RADIUS 


Method I 


The data presented in the preceding section 
showed that the size of the sedimentation radius 
departs from that of the true nominal radius 
with increasing Reynolds number. For many 
practical purposes it is, however, of interest to 
obtain a size classification of granular substance 
by sedimentation or elutriation in air or water, 
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although the computed size values depart from 
the actual ones with increasing R-values. Only 
when dealing with perfect spheres is the true 
volume size obtained, because then the actual 
radius, 7., equals 7, and 7,. Stokes’ or Oseen’s 
well-known formulas, the validity of the latter 
extending to an R-value of about 0.5, are used 
for lower Reynolds numbers. A recent improve- 
ment by Goldstein” carries the validity of his 
formula but slightly above Reynolds number 0.5. 
Zahm’s'® equation for the coefficient of resistance 
for spheres at higher Reynolds numbers is valid 
between R-values ranging between approxi- 
mately 1000 and 0.3, but, instead of gradually 
merging with the linear function of Stokes’ law, 
the curve constructed on his formula departs 
from Stokes’ law in an increasing degree with 
decreasing Reynolds. This deviation, as well as 
the fact that the radius must be known in order 
to compute the value of the coefficient of re- 
sistance, makes Zahm’s formula unsuitable for 
sedimentation purpose. There are several em- 
pirical formulas which hold approximately for 
higher Reynolds numbers, but they are usually 
not confined to any special shape and disregard 
the influence of the form under the author’s 
assumption that the shape variation of sedi- 
mentary particles is slight, or else they fail to 
bridge satisfactorily the region between the upper 
limit for the validity of Stokes’ law and higher 
Reynolds numbers. 

Every factor in Eqs. (2) and (3) must be 
known in order to determine the Cre=/(R) in 
the logarithmic chart. The size of the sedimenta- 
tion radius can, however, be determined without 
knowing the size of the true nominal radius. 
An inspection of Eqs. (2) and (3) discloses that 
the true nominal radius, 7, enters both formulas. 
By combining the two into a ratio, 0, the true 
nominal radius, 7,, is eliminated : 


0=Cr/R=4(ps— pr)gv/3v* pr. (4) 


We assume that a disk-shaped particle No. 1 
and a spherical particle No. 2 have the same 


17S, Goldstein, The Steady Flow of Viscous Fluid Past a 
Fixed Spherical Obstacle at Small Reynolds Numbers, 
Proc. Roy. Soc. (London) A123, 225 (1929). 

1A. F. Zahm, Flow and Drag Formulas for Single 
Quadrics, Twelfth Annual Report of the National Advisory 
Committee for Aeronautics, 1926, Report No. 253, pp. 517- 
537 
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specific gravity and the same terminal uniform 
settling velocity in the same fluid. Thus all 
properties of the two sedimentation systems are 
the same, except the true nominal radii of the 
two solids, which we can suspect to differ on 
account of the same velocity of the two solids, 
but the relative great difference in shape between 
a sphere and a disk. Since the true nominal 
radius, f,, is eliminated in the ratio o, and 
because all the other properties of the two 
systems are the same, we conclude that the ratio 
0, of the disk equals 02 of the sphere. Since the 
velocities of the two solids were assumed to be 
the same, we can formulate the following thesis: 
Thesis I. All solids of the same specific gravity, 
irrespective of their volume sizes, shapes or shape 
positions, have the same terminal uniform settling 
velocity in the same fluid, when the ratio of the 
coefficient of resistance to the Reynolds number for 
respective solids has the same numerical value. 
The sedimentation radius can be determined 
graphically on the basis of this thesis. All solids, 
which according to their Cg =f(R) coincide with 
the linear function y=x or Ceg=R, in the special 
case of chart Fig. 4, must have the same terminal 
settling velocity under the provisions of Thesis I, 
because the ratios y/x or Cr/R have the same 
numerical value, i.e., one. The same relationship 
exists between solids, the Cr=f(R) of which 
coincide with any other straight line parallel to 
y=x. Consequently, disk No. 1 (chart Fig. 4) 
on line a to b||y=x has the same terminal 
velocity and the same sedimentation radius as 
a sphere No. 2 on the same line and on the 
experimentally constructed curve for spheres. 
The data required for our numerical illustra- 
tion have been obtained from Schmiedel.’ Dot 
No. 1 in chart Fig. 4 represents the Cr=/(R) 
for a disk of steel (y¥=0.1219) in its tardiest 
shape position. The coefficient of resistance and 
the Reynolds number have been computed 
according to Eqs. (2) and (3). The numerical 
values are given in Table I, solid No. 1 
(Schmiedel 13). In order to construct a case in 
analogy with ordinary sedimentation methods, 
we assume that the Cr=/f(R) and the true 
nominal radius of the disk are not known. We 


‘compute the ratio 0; for the disk by applying 


Eq. (4) and substituting the numerical values of 
Table I. 
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4(7.674 —1.199)981 -0.3684 
0);>= = 
3-(11.55)3- 1.199 





1.6889. (5) 


Table I shows that Cr = 26.393 and R= 15.627, 
“. Cr/R=26.393/15.627 =1.6889, which is the 
same numerical value obtained by Eq. (5). 
The numerical value of 0; does not permit any 
determination of the Cr=f(R), but we know 
that, if the numerical value of o for any other 
solid equals that of 0; under the provision of 
Thesis I, then this solid has the same terminal 
velocity as the disk. In order to construct the 
line a to b we choose an arbitrary coefficient of 
resistance, for instance, Ce =40, and determine 
the Cre=f(R) and the Reynolds number for a 
point a of the same numerical o-value as the 
disk. Thus 40/R=1.6889, R=23.68. The 
cross-mark a in chart Fig. 4 is located by 
Cr=40 and R= 23.68. We can now either draw 
a line || y=x, or locate a point b in the same 
manner as point a was obtained. The line a to b 
crosses the experimentally constructed curve for 
spheres at a point No. 2, which is the Ce =f(R) 
for a sphere of the same terminal velocity as the 
disk, i.e., the sedimentation sphere of the disk. 
Reynolds number 4.50 for point No. 2 is de- 
termined at the lower margin (x-axis) of the 
chart. The sedimentation radius, r,, can now be 
computed by Eq. (3). Substituting the numerical 
values of Table I and inserting r, instead of r, 
in Eq. (3), we find that: 


4.50 = (2-7r,- 11.55) /0.3684, 
.. £2 =0.071767 cm. 


(6) 


The accuracy of the obtained value of 7, 
depends on the exactness of the reading of the 
Reynolds number on the chart. However, even a 
crude determination of the R-value is satis- 
factory, because the variation of the experi- 
mental data, upon which the curve for spheres is 
based, permits a sufficient range in variation of 
the determined Reynolds number (see Fig. 5). 


Method II 


The well-known equation for the Stokes radius, 
7s:, and the Stokes velocity, vs1, read: 


rse=(9uv/2(ps—pr)g ]}. 


Vse=2r(ps— pr)g/9u. 


(7) 
(8) 
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For spherical shapes within the region of 
Stokes’ law the actual radius, r., equals the 
Stokes radius, rs:;, and the Stokes velocity, vs,, 
equals the actual velocity, va. For non-spherical 
solids in the tardiest position the true nominal 
radius, fn, is always larger than the Stokes 
sedimentation radius, 7,s:, computed according 
to Eq. (7). For R-values higher than 0.2, rs:<ra 
for spherical shapes, and for non-spherical forms 
in the tardiest position r.5:<7n, and %,<?fn, 1s 
representing the sedimentation radius computed 
according to Method I. 

The linear function of Stokes’ Jaw is expressed 
by the equation Crs,=24/Rs:, where Crs: 
represents the Stokes coefficient of resistance and 
Rs, is the Stokes Reynolds number. The formula 
is obtained by inserting Stokes’ well-known 
expression for the resistance in the conventional 
formula for Cr: 














R OrursWse 
Cr= ; Cesr= < ’ 
mr?(v* pr) /2 trsi(vsrpr)/2 
12u 
RSst= (9) 
’stUstPr 

Rs:= 2rssipr/p= 2rsWst/v. (10) 
CresiRsi= 24, Crst= 24/Rs:. (11) 


Since rs:<7qand vs; >. for Reynolds numbers 
higher than 0.2, Eq. (11) cannot, strictly speak- 
ing, be employed for computing the Stokes 
coefficient of resistance for higher R-values, 
because the numerical values of r. and v, have 
not the same relation to each other as 7s, and 
vst, the latter being in agreement with Stokes’ 
law, according to which the square of the radi 
is proportional to the velocities. Eq. (11), or 
rather a similar one, Cge=24/R, has, however, 
frequently been used for computing the Stokes 
coefficient of resistance for higher R-values, with 
the purpose of showing the amount of deviation 
of the Stokes coefficient of resistance from the 
actual one, Cra. Since rgsi:<fa, Use>Va, Cree 
<Cra, for Reynolds numbers higher than 0.2, 
we have really two different conceptions in- 
volved, namely, the Crsi=f(Rsm) correspond- 
ing to Cra=f(R.) on the basis of vg,=va, and 
the Crsz=f(Rsi2), which corresponds to Cra 
=/(R.) on the basis that rg;=7a. 





f 
: 
| 
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*. Use=Va, ..» Crs =12u/rsWapr. (12) 
- Use=Va, «|. Rou =2rsWapr/u 

=2rsWa/v. (13) 

Crsu=24/Rsu. (14) 

“Prst=fa, -» Crsie=12pu/ratsipr. (15) 
“Prst=fa «. Rsgw=2radsipr/u 


=2rvs./v. (16) 
Crsi2=24/Rsv2. (17) 


Since va<vs; and ra>rg:, hence Cra>Crsi 
>Crse and Rsgin<Ra<Rsgiz. The conceptions 
are graphically demonstrated in chart Fig. 5, 
where dot No. 3 represents the Cra=f(R.) fora 
sphere, and dots Nos. 1 and 2 represent Crsi1 
=f(Rsu) and Crswz=f(Rsi2), respectively, for 
the same sphere in the same fluid. Our aim is to 
show that the actual radius, r,, and the actual 
sedimentation radius, 7,, can be determined 
graphically by correction of the Stokes radius, 
rs:, and the Stokes sedimentation radius, r,s,. 
We have chosen a numerical illustration as the 
simplest manner of demonstrating the method. 

Applying Eq. (7) and substituting the nu- 
merical values of Table I, solid No. 1 (Schmiedel 
13), we find that the Stokes sedimentation 
radius, 7,s:, for the disk amounts to, 


9-0.4417-11.55 
rest= =0.06012 cm. (18) 
2(7.674—1.199)981 





Next we determine the Reynolds number, 
Rs, for the Stokes sedimentation radius by 
formula No. 13, and svbstituting numerical 
values: 


Rosi = (2-0.06012- 11.55) /0.3684=3.7697. (19) 


We can now compute the Stokes coefficient of 
resistance, Crs, by Eq. (14). 


Crsu=24/3.7697 = 6.367. (20) 


Dot No. 3 in chart Fig. 4 has a Ces: = 6.367 
and a Rs: =3.7697, thus located on the point 
where the linear function of Stokes’ law crosses 
the line a to b, the latter crossing the experi- 
mentally constructed curve at a point No. 2, 
i.e., the point which by Method I was determined 
as the Cra=f(R.) for the disk’s sedimentation 
sphere of radius r,. The construction of the line 
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a to b and subsequent computations are the same 
as in Method I. Thus we find that r,=0.07176 
cm, which is the same numerical value obtained 
by Method I. 


COMPUTATION OF THE ACTUAL RADIUS OF A 
PERFECT SPHERE AND OF THE SEDI- 
MENTATION RADIUS 


Method III 


We have previously stated that the curve for 
spheres in Figs. 1, 2, 3 and 4 was constructed on 
the basis of experimental data, corrected by 
Ladenburg’s formula for the influence of the 
wall of the sedimentation vessel. The value of 
the computed sedimentation radius refers there- 











TABLE II. 
I II Ill IV 
Rsn Crsu CRa R, 
0.010 2400.00 2407.68 0.01003 
0.100 240.00 243.84 0.1016 
0.250 96.00 98.91 0.2575 
0.300 80.00 82.75 0.3103 
0.350 68.57 71.20 0.3634 
0.400 60.00 62.52 0.4168 
0.500 48.00 50.36 0.5246 
0.600 40.00 42.29 0.6335 
0.700 34.28 36.42 0.7436 
0.800 30.00 32.05 0.8547 
0.900 26.66 28.64 0.9668 
1.000 24.00 25.92 1.080 
1.500 16.00 17.69 1.659 
2.000 12.00 13.55 2.259 
2.500 9.600 11.05 2.879 
3.000 8.000 9.379 3.517 
4.000 6.000 7.264 4.843 
5.000 4.800 5.982 6.232 
6.000 4.000 5.119 7.679 
7.000 3.428 4.497 9.182 
8.000 3.000 4.026 10.73 
9.000 2.666 3.657 12.34 
10.00 2.400 3.360 14.00 
15.00 1.600 2.449 22.96 
20.00 1.200 1.979 32.98 
25.00 0.9600 1.688 | 43.97 
30.00 0.8000 1.489 55.86 
40.00 0.6000 1.232 82.16 
50.00 0.4800 1.071 111.6 
60.00 0.4000 0.9597 143.9 
70.00 0.3428 0.8772 179.1 
80.00 0.3000 0.8133 216.8 
90.00 0.2666 0.7621 257.2 
100.0 0.2400 0.7200 300.0 
150.0 0.1600 0.5848 548.2 
200.0 0.1200 0.5096 849.3 
250.0 0.09600 0.4602 1198.6 
300.0 0.08000 0.4248 1593.1 
400.0 0.06000 0.3762 2508.2 
500.0 0.04800 0.3436 3580.0 
600.0 0.04000 0.3199 4798.3 
700.0 0.03428 0.3014 6155.2 
800.0 0.03000 0.2866 7644.2 
900.0 0.02666 0.2743 9260.8 


























SOME NEW 


fore to the corrected values for spheres. We shall 
in the following present a formula, by which the 
sedimentation radius and the actual radius of a 
perfect sphere can be computed directly, without 
graphic methods and on the basis of experi- 
mental data without Ladenburg’s correction. 

In chart Fig. 5 the dotted line, J to 3, con- 
necting the linear expression for Stokes’ law 
with the experimentally constructed curve for 
spheres of mnon-corrected data, represents the 
ratio o of a given numerical value (0 =0.0914). 
Thus a sphere of a radius rs (point No. 1, 
Fig. 5) must have the same terminal uniform 
settling velocity’? (vgsi::=v.a3) in the same fluid 
as a sphere of the same specific gravity and of an 
actual radius, 743; the Cra=f(R.) of the latter 
sphere being located on the experimentally con- 
structed curve at a point No. 3 of the same o- 
value as point No. 1. Since all properties of the 
“two sedimentation systems,” except the radii, 
are the same, i.e., PS1>PS3, PF1>—PF3, Mi> Hs, 
£1= 23, Vsti =Vaz, it follows that, 


Cras 8fas(pss— prs) £3 

















Crsua 3va3" PFs 
30 su’ pri Yas 
amon, (21) 
8rsu(psi- pri) £1 st 
Ras 2rasVas VY} Ta3 
= ; Bam, (22) 
Rsn V3 2rsutsa Tsu 
Ra3/Rsu= Cra3 /Crsu. (23) 


It is possible to carry these equations a step 
further by determining the numerical constant 
of either Ra/Rsin or Cra/Crsi. We have found 
by statistical methods that, 


R./Rs11 =1+0.08(R 543) 9897. (24) 
R.=Rsu[1 +0.08(R 11) ©9897 J. (25) 
Cra= Crsul1+0.08(Rs5 11) 597], (26) 


where the exponent 0.69897 =log 5. 
Column I of Table II gives successive numbers 


19 The reasoning above is based on the difference in the 
graphic positions of the two spheres, a difference due to the 
failure of Stokes’ law for higher Reynolds numbers. The 
reasoning is not obscured by the fact Cr=f(R) of all 
spheres really should coincide with the experimentally 
constructed curve. 
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of Rsu, and column II gives corresponding 
Crsu according to Stokes’ law, Cres =24/Rsu. 
Columns III and IV give the actual coefficient 
of resistance, Cra, and the actual Reynolds 
number, R,, respectively, for Crsu and Rgz. 
The values of R, and Cr, have been computed 
according to Eqs. (25) and (26), respectively. 

The solid curve in Fig. 5 illustrates the com- 
puted values of Cra and Ra, expressed as Cra 
=f(R.). The experimental data cluster round the 
curve with sufficiently close coincidence to 
verify the formulas. The curve merges with 
those of Oseen’s and Goldstein’s formulas at an 
R-value of about 0.5. It approaches the linear 
function of Stokes’ law with decreasing Reynolds 
numbers, and, although the curve theoretically 
never crosses the Stokes line, it practically 
merges with the linear function of Stokes’ law 
at an R-value of 0.01 (see Table II). The curve 
appears to depart somewhat from the experi- 
mental data at an R-value of about 3500. At 
these high Reynolds numbers ordinary sedimenta- 
tion apparatus are too small for measuring the 
velocity in water. Ordinary sedimentary particles 
have, at this high R-value for water as sedi- 
mentation fluid, sufficient size to permit separate 
weight determination by ordinary laboratory 
balances. If, however, sedimentation methods 
are used for R-values above 10,000, the sedi- 
mentation radius can be computed by Method I. 

The actual radius, ra, of a sphere and the 
sedimentation radius, 7,, for a non-spherical 
solid can be computed in various ways. We have 
chosen the simplest method. Eqs. (13), (22) 
and (24) show that, 


R./Rsu = ra/t sti =1 +0.08(R gs 41) ©9897, 
Hence, 


rs=la="siL1+0.08(27s wapr/u)*** ], 


where 


(27) 


's¢= [9uv./2(ps ans pr)g |}, 


v2 is the observed, terminal uniform settling 
velocity, 7. the actual radius of a perfect sphere, 
and r, the sedimentation radius for a non- 
spherical solid. 

Table III gives the computed and measured 
radii of spheres for a range of Reynolds numbers 
down to an R-value of 0.0531, the latter being 
well below the upper limit for the assumed 
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TABLE III. 

Observed Reynolds Measured Computed (Stokes) Error Error 
Data by Ps Pr mm Va number Ve Fs rst % q 

cm/sec. R, cm cm cm (Wadell) (Stokes) 
Allen 7.666 0.99963 0.012735 100.5 3144.5 0.1993 0.2042 0.029676 +2.4 —85.1 
Schmiedel 7.818 1.227 1.3718 36.24 17.922 0.27645 0.27081 0.18601 —2.03 —32.7 
Schmiedel 2.700 1.228 1.6246 1.666 0.2437 0.09680 0.09448 0.091840 -—2.39 — 5.12 
Schmiedel 2.700 1.233 1.9987 1.370 0.1636 0.09680 0.09455 0.092530 -—2.32 — 4.41 
Schmiedel 7.524 1.233 2.0452 1.703 0.1059 0.05045 0.05122 0.05040 +1.52 — 0.09 
Schmiedel 2.700 1.228 1.5939 0.6109 0.0531 0.05655 0.05564 0.055085 -—1.60 — 2.59 











validity of Stokes’ law. The data, illustrating the 
general difference between the results obtained 
by Eq. (27) and by use of Stokes’ law, refer to 
spheres, the Cp =f(R) of which coincide as close 
as possible with the constructed curve. The 
difference between the computed and measured 
radii is expressed as an “error,’’ computed as 
percentage of the value of the measured radius. 
The plus sign (+) shows that the computed 
radius is larger than the measured one, and the 
minus sign (—) indicates that the former is 
smaller than the latter. The table shows that 
even in the vicinity of and below Reynolds 
number 0.2 the formula presented in this paper 
gives results, which agree well with the experi- 
mental data. 


A NEw RESISTANCE FORMULA FOR SPHERES 


For convenience we repeat the conventional 
resistance formula for spheres: 


R = Tra? (Va2 pr/2)Cra- (28) 


Eliminating Cr. and substituting its equiva- 
lent given in Eq. (26): 
R=nr2(v.2pr/2) 
x Cersull +0.08(R 511) 97 J. 
Eliminating Ces and substituting its equiva- 
lent given in Eq. (12): 
R = (rr.0a?pr/2) si (1 2u/rs apr) 
x [1+0.08(27r5 wapr/m) 97 J. 


Hence the following equation for the resist- 
ance, R; 


(29) 


(30) 


R=6rprata: (%a/1 st) 
<[1+0.08(27s vapr/u) 7], (31) 
where v, is the observed terminal velocity, and 


fs, is computed in the same manner as in 


Eq. (27). 


COMPUTATION OF THE TERMINAL UNIFORM SET- 
TLING VELOCITY OF A SPHERE OF A RApDIUS 
EQUAL TO THE TRUE NOMINAL RADIUS 
OF A GIVEN NON-SPHERICAL SOLID 


By combining Cr and R (Eqs. (2) and (3)) 
into a product, p= Cr(R)’, the velocity factor, v, 
is eliminated : 


b= Cr(R)* = 327,3(ps— pr)g/3prv’. 


We assume that a disk-shaped particle No. 1 
and a spherical particle No. 4 have the same 
true nominal radius and the same specific gravity, 
and that they are settling in the same fluid. 
Thus all properties of the two sedimentation 
systems are the same, except that the shapes of 
the solids differ and the terminal velocities are 
not known. Applying Eq. (32) we obtain two 
products, p; and f for the disk and the sphere, 
respectively. Since the velocity factor is elimi- 
nated, and since all corresponding factors of p; 
and #, equal each other, we conclude that the 
product p; for the disk equals ~, for the sphere. 
We can then formulate the following thesis: 

Thesis II: All solids of the same specific gravity 
settling in the same fluid, irrespective of their 
shapes, shape positions and terminal settling 
velocities, have the same true nominal radius or 
the same volume size, when the product Cr(R)* of 
respective solids has the same numerical value. 

In chart Fig. 4 we have a straight line repre- 
senting the equation y= 100/x?, or in the special 
case of this chart, Cre=100/R?, hence, Cr(R)’ 
= 100. Consequently, all solids the Ce=f(R) of 
which coincide with the line y=100/x? must 
have the same numerical p-value, i.e., p= Cr(R)* 
=100, and they must therefore also have the 
same true nominal radius, 7,, under the pro- 
visions given in Thesis II. The same relationship 
exists between solids, the Ce=f(R) of which 


(32) 
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coincide with any other straight line || y 
= 100/x2. 

We have chosen a numerical illustration as 
the simplest method of demonstrating the com- 
putation of the terminal velocity for a sphere of 
a radius equal to the true nominal radius of a 
given non-spherical solid. Dot No. 1 in chart 
Fig. 4 represents the Ce=f(R) for a circular 
disk of steel in its tardiest shape position. The 
known properties of the sedimentation system 
are, fn, ps, pr, and v. Applying Eq. (32) and 
substituting the numerical values of Table I, 


solid No. 1 (Schmiedel 13): 
32-(0.24922)8- (7.674 —1.199)981 





P\= = 6445. (33) 


3-1.199- (0.3684)? 


If the numerical p-value of any other solid 
equals p; under the provisions given in Thesis II, 
then this solid must have the same true nominal 
radius as the disk. We choose an arbitrary 
coefficient of resistance, for instance, Cr =40, 
and then determine the Reynolds number for 
this’ solid h: 


pi: =6445=40(R)?, .. R=12.69. (34) 


The cross-mark h in chart Fig. 4 represents 
the solid h of Cp =40, and R=12.69. This solid 
must have the same true nominal radius as a 
sphere of a radius equal to the true nominal 


radius of the disk. It follows also that the 
Cr=f(R) of the disk and of the sphere must 
coincide with a straight line passing through 
h| | =100/x*. In construction of such a line we 
choose an arbitrary coefficient of resistance, for 
instance, Cre=3, and determine the point gq in 
the same manner as / was located. We extend 
the line h to q till it crosses the experimentally 
constructed curve for spheres at a point No. 4. 
The Reynolds number for this point may be 
read at the lower margin (x-axis) of the chart. 
Having found that the Reynolds number for 
sphere No. 4 amounts to 72, we can compute 
the terminal velocity for this sphere. Applying 


Eq. (3) and substituting the numerical values of 
Table I: 


R=72=2-0.24922-0,/0.3684, 


35 
‘ .. va =53.216 cm/sec. (35) 


We have thus arrived to the result that a 
sphere of the same volume size and of the same 
specific gravity as the disk has a terminal 
uniform settling velocity of 53.216 cm/sec. in 
the same fluid. This velocity is about 4.6 times 
greater than that of the disk (11.55 cm/sec.) in 
its tardiest shape position. 

The actual velocity can also be computed by 
correction of the Stokes velocity on the basis 
of the same reasoning presented above. 
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Thermal Agitation Voltages in Resistors 
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The limited range of previous measurements of thermal agitation voltage is pointed out 
and data are given which extend the range of resistances and frequencies considerably. The 
results substantiate Nyquist’s theory for resistances up to 1.7 X 10° ohms and for frequencies 
between 1 and 10,000 cycles per second. A value of 1.366 X 10-" ergs per degree with a probable 
error of 0.016 10~" ergs per degree is obtained for Boltzmann's constant. This is nearer to 
the accepted value than previous determinations by this method have given. 





OR some years it has been known that the 

thermal motion of electrons in conductors 
gives rise to a voltage across the terminals of 
the conductor. The effect was discovered and 
first studied by J. B. Johnson.' This voltage 
may, in the case of a high-gain amplifier, be 
sufficiently large to mask the voltage to be 
measured and therefore determine the lower limit 
of measurements. It is the purpose of this paper 
to discuss the effect and its relation to the 
problem of high-gain amplification and to present 
additional data as to its nature. 

The nature of the phenomenon may be easily 
explained in very simple words. In every con- 
ductor there exists free electrons. These elec- 
trons, because of their thermal energy, are in a 
continual state of agitation. The movement of 
these electrons constitutes a current, and when 
a current flows through a resistance a difference 
of potential exists across the resistance. 

Certain characteristics may be immediately 
pointed out upon the recognition of the existence 
of the effect. First, the voltage is a random 
voltage since it is caused by the random motion 
of the electrons. Second, the average value of 
the voltage over a sufficiently long period of 
time is zero because there is no general trend to 
the motion of the electrons, that is, there is no 
piling up of the charge at any point in the 
conductor. Third, the voltage is a function of 
the temperature and the resistance. 

By assuming that the ordinary laws of con- 
duction hold, it is possible to obtain additional 
data on this effect. Consider two resistances R; 
and R, of any magnitude, physical size, shape or 
composition connected as shown in Fig. 1. Let 
é: be the mean-square voltage generated in 


1 Johnson, Phys. Rev. 32, 97 (1928). 


resistor R, and e,2 that generated in Re. By 
virtue of the voltage generated in R; a current 


iy =e:2/(Ri+Rs) (1) 


will flow in the circuit. The rate at which energy 
from R;, is dissipated in R, is 


P, =i;?Re =eRe/(Rit+R2)®. (2) 


Similarly, the rate at which energy from R; is 
dissipated in R, is 


P.=ie?R, = e22R,/(Ri+R2)?. (3) 


But, in order that there be no net energy 
transfer, that is, that one resistor shall not be 
cooled and the other heated, these powers must 
be equal and e,?R2= eR, or 


er, e2= R; ‘Ro. (4) 


We therefore conclude that, if the ordinary laws 
of conduction hold, the mean-square value of 
the voltage is directly proportional to the re- 
sistance and is independent of the physical size, 
shape, and composition of the resistor. H. 
Nyquist,” by replacing R, of Fig. 1 by a complex 
impedance R;+jX, and providing that there be 
no net energy transfer, has shown that for a 
complex impedance the mean-square voltage is 
proportional to the real part of the impedance. 
In amplifiers the principal circuit encountered 
is that of a resistance and capacitance in parallel, 
the capacitance being that of the circuit leads 
and of the tube. For such a circuit the real 
part of the impedance is 


R(f) =R/(1+R°C*e"). (5) 


2 Nyquist, Phys. Rev. 32, 110 (1928). 
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Fic. 1. Circuit for showing the variation of resistance noise 
with resistance. 


J. B. Johnson’s! work indicated that the mean- 
square voltage was proportional to the tempera- 
ture and the resistance and independent of the 
frequency except insofar as the frequency affected 
the resistance. Since an amplifier is required to 
measure the voltage it becomes convenient to 
express the effect in terms of the output current 
in an amplifier. Johnson found that if Y(f) is 
the conductance function of the amplifier the 
output current is expressed by 


Paer f “R(F)| YP) Pf, (6) 


where k is Boltzmann’s constant and T is the 
absolute temperature. His measurements covered 
temperature from — 180°C to 100°C, frequencies 
from 300 to 2000 cycles per second, and re- 
sistances up to almost one megohm. Using Eq. 
(6) he obtained a value of 1.27 10~"* ergs per 
degree for Boltzmann’s constant. H. Nyquist? 
developed a theoretical formula and also arrived 
at Eq. (6). 

H. D. Ellis and E. B. Moullin® used the effect 
as a means of determining Boltzmann’s constant, 
assuming Eq. (6) to be correct. Their measure- 
ments were principally limited to two values of 
resistance, 62,000 and 110,000 ohms, and covered 
the frequency range between 1700 and 3300 
cycles per second. They give values of Boltz- 
mann’s constant ranging from 1.35X10-* to 
1.361 10-" ergs per degree. 

E. K. Sandeman and L. H. Bedford‘ made 
measurements with resistances between 104 and 
10° ohms but for wide-frequency band. In one 
case they used a band between 30 and 7000 
cycles per second and in another case a band 
between 30 and 2200 cycles per second. 

Williams and Wilbur’ reports measurements 

§ Ellis and Moullin, Proc. Phil. Soc. 28, 386 (1932). 


*Sandeman and Bedford, Phil. Mag. 7, 774 (1929). 
5 Williams and Wilbur, Science 76, 519 (1932). 


VOLTAGES 293 
made at boiling water and liquid air temperatures 
but are not clear as to the range of resistance 
and frequency covered. 

Although the above researches give good 
checks on the theory, their range is somewhat 
limited. The use of comparatively small re- 
sistances and low frequencies made it unneces- 
sary to take into account the variation of re- 
sistance with frequency. In addition, work 
reported by Thatcher® seemed to indicate that 
the theory did not hold for higher resistances. 
It was also thought possible that an abnormal 
effect might occur at very low frequencies as in 
the case of the small shot effect. 

Measurements were therefore made for re- 
sistances between 1 X10* and 1.7 X 10° ohms and 
at frequencies between 1 and 10,000 cycles 
per second. The measuring circuit was arranged 
to make the shunt capacitance as smal] as 
possible. This was desirable since the size of the 
capacitance determines the range of resistance 
that may be measured at a given frequency. 
Curve I of Fig. 2 gives R(f)/R=1/(1+ R?C*w*) 
as a function of RCw. For values of RCw less 
than 0.1, R(/) is equal to R to within 1 percent. 
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Fic. 2. Variation of the resistance of a parallel R—C 
circuit with frequency and d.c. resistance. 


6 Thatcher, Phys. Rev. 40, 121 (1932). 
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Fic. 3. Arrangement of the first stage and resistor mounting. 
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Fic. 4. Circuit diagram of the first stage and calibration 
circuit. 


However, for values of RCw greater than 0.1, 
R(f) may be considerably different from R. 
Curve II of Fig. 2 gives R(f{)Cw as a function 
of RCw or R(/) as a function of R for a fixed 
capacitance and frequency. It reaches its maxi- 
mum value when RCw=1 or when R=1/Ca, 
and at this point R(f) =0.5 1/Cw. The maximum 
resistance obtainable at a given frequency is 
therefore determined by the input capacitance 
of the amplifier. Fig. 3 shows the physical 
arrangement of the first stage. The input capaci- 
tance as measured on an impedance bridge was 
13 micromicrofarads. 

The amplifier used was an ordinary resistance- 
capacitance coupled amplifier having a tuned 
circuit in the plate circuit of one stage. The 
first stage was located in a separate shielded box 
with the measuring circuit. Fig. 4 is a schematic 
diagram of the first stage showing the arrange- 
ment of the resistor upon which measurements 
were made and of the calibration circuit. . 

Errors in the calibration circuit due to shunt 
capacitance were eliminated by shielding. Fig. 5 
is a schematic diagram of the calibration circuit 
showing the shunt capacitances and the im- 
pedance R, of the cable from the voltmeter. 

Condenser C; causes an error only because of 
the additional current it draws through R,. 
Since R, is only a small fraction of an ohm, this 
error is negligible. Condenser C2 is in parallel 
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Fic. 5. True diagram of the calibration circuit. 


with 1 ohm and would have to be two micro- 
farads in order to cause an error of 1 percent at 
10,000 cycles per second. Condenser C; is in 
parallel with 1 megohm and is apt to cause an 
error at high frequencies unless it is made small. 
By shielding R; from Re, C; may be reduced to 
the self-capacitance of R;. Wire-wound resistors 
having a self-capacitance of about two micro- 
microfarads may be used with an error of less 
than 1 percent for frequencies up to 10,000 
cycles per second. By shielding R; and the cable 
to the amplifier, errors due to inductive coupling 
are eliminated. 

The output of the amplifier was fed into a 
thermal meter so that the root-mean-square 
value of the output current could be read. 

Johnson’s data were considered sufficient evi- 
dence of the relation of the voltage to tempera- 
ture so that measurements were made only at 
300° absolute with a probable error of 1°. 

Measurements were carried out on the as- 
sumption that the present theory did not hold 
so that precautions were taken to eliminate 
reliance on the formula of Eq. (6). It was 
assumed that the voltage was a function of only 
the resistance, the temperature and the fre- 
quency. Then in the small increment of frequency 
between f and f+d/ the current in the output 
of the amplifier would be given by 


di?= F(R(f), T, f]| Y(f) | df, (7) 


and the total current is given by 


ia f F(R(f), T.f)| YP) df. (8) 


Now if the amplifier is tuned very sharply, F 
may be considered constant over the frequency 
range that contributes appreciably to 7 and 
Eq. (8) becomes 


i2= F(Ro, T, fo) f "| Y(f) |*df, (9) 














where Ro and fo are the values of R(f) and f at 
which Y(f) is a maximum. The band width 
varied from about 10 cycles at 10,000 cycles per 
second to about 1 cycle at 30 cycles per second 
as compared with band widths ranging between 
200 and 400 cycles used by Ellis and Moullin. 

Now on the assumption that 7 is proportional 
to T Eq. (9) becomes 


P=GRAIT{ \YU)Idf. (10) 
“se 

The method of measurement consisted in 
tuning the amplifier to a given frequency and 
then reading the current for a series of resistances 
up to 250 megohms, and also for zero resistance. 
The conductance function, | Y(f)|, was then 
obtained by applying a variable frequency known 
voltage to the calibration circuit. After correction 
is made for the noise from the amplifier, the first 
of these readings gives a series of values of 7? in 
Eq. (10) for various values of Ro. The second 
gives the value of the integral* of Eq. (10). 
From these results we can obtain the value of 
G as a function of Ro for a given frequency. 

Fig. 6 shows (7)!, of Eq. (10), typically 
plotted against R and Fig. 7 gives the corre- 
sponding curves of | Y(f)| and | Y(f)|*% The 
narrow band width was made possible by con- 
necting a dynatron in parallel with the tuned 
circuit. 

The slope of the curve of Fig. 6 is equal to 2 
on the straight part of the curve, indicating that 
G varies directly as Ro, which would be indicated 
by Eq. (6). If the curved portion is to be ac- 
counted for by shunt capacitance, then the 
value of the capacitance may be calculated from 
the frequency and the resistance at which (7)! 
becomes a maximum. At this point R=1/Cw or 
C=1/Rw. This gives C=12.7 micromicrofarads 
as compared to 13 micromicrofarads, the meas- 
ured value. The solid line of Fig. 6 gives the 
theoretical value assuming that 7? varies as the 
real part of the impedance, that is, G of Eq. (10) 
is directly proportional to Ro. The only point 
that is off an appreciable amount is the one at 
2.5 X 108 ohms and it is probably an experimental 
error. 

Measurements were made at frequencies of 
approximately 30, 50, 100, 160, 250, 500, 780, 


* See appendix. 
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Fic. 7. The conductance and conductance squared of the 
amplifier for one frequency used. 








0 
3170 


1000, 2000, 3000, 4000, 7000 and 10,000 cycles 
per second. 

In addition to the above frequencies measure- 
ments were made at approximately 1 cycle 
per second and with resistances as high as 1.7 
10° ohms using resistance capacitance tuning 
throughout the amplifier. The results at this 
frequency are not accurate because it was 
necessary to use a rectifying meter instead of a 
thermal meter in the output of the amplifier. 
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TABLE [. 
(4%)? —m.a. Appar. 
Freq Ps VY (f) *df —(m.a./V)? for input 
for Calcu- Meas- Total Ro = 108 capac- 
max. Y lated ured ohms K X 10% itance 
1. 2. ae 4. 5. 6. A 
0.85 0 206 206 1.75 1.24 —_ 
29.32 126 941 1067 4.24 1.41 —_ 
51.53 26 467 493 2.79 1.31 — 
100.66 80 1288 1368 4.74 1.37 13.3 
159.46 590 7570 8160 11.2 1.28 13.3 
250.66 4200 42000 46200 28.4 1.46 11.5 
500.7 8200 87000 95200 40.3 1.42 12.7 
773.3 8700 62300 71100 35.2 1.45 12.1 
998.0 35000 212000 247000 60.4 1.23 12.3 
2042.5 6500 635000 70000 34.8 1.44 13.0 
3139.6 8900 90000 98900 40.4 1.38 12.7 
3994.0 4200 100600 104800 41.0 1.34 13.5 
7173.5 1100 52800 53900 26.3 1.07 13.0 
10640.0 2800 45600 48400 28.0 1.35 12.5 
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However, the measurement is of sufficient im- 
portance to include with the other data as it 
indicates no change in the theory at very low 
frequencies. 

Table I is a tabulation of the results at the 
various frequencies. Column 6 gives the value of 
Boltzmann’s constant as computed from Eq. 
(6). By neglecting the values for 1 and 7000 
cycles per second, these data give an average 
value of 1.366X10-" erg per degree with a 
probable error of 0.01610-"* erg per degree. 
Fig. 8 is a curve of k plotted against frequency. 
The points are scattered at random about the 
accepted value with no definite trend. 

Column 7 of Table I gives the apparent input 
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Fic. 8. The measured value of Boltzmann's constant as a 
function of frequency. 
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Fic. 9. Equivalent circuit of the stage containing the 
tuned circuit. 





capacitance as computed from the value of d.c. 
resistance and the frequency for maximum (7)!. 

These data seem to indicate that the present 
theory is true for the range of resistance and 
frequency covered in this research. 


APPENDIX 


Since the total area under the | Y|?—f curve would 
have been very difficult to measure, the portion far re- 
moved from the resonant frequency was calculated by an 
approximate method. 

In Fig. 9, R, represents the plate resistance of the tube 
and Z the impedance of the tuned circuit. The voltage 
amplification of this circuit is 


e2/e,=u!Z!|/|Rp+Z). (11) 


If R, is large, then sufficiently far from resonance |Z|<R, 
so that 


e:/e2="|Z|/Rp (12) 


or is directly proportional to |Z|. Neglecting the error 
caused by the frequency characteristic of the remainder 
of the amplifier, it is true that | Y| is also proportional to 
|Z|. By setting 


|Z| =Lw/(1—LCw*) =(Bflfo)/L1—(f/fo®] = (13) 
we obtain 
| ¥| =B’x(f/fo)/L1—C/fo)*] 
and 
Boas > (B’f/fo)* ‘ ada 
|} \°df = | ——————_df = Bf. |} ————_,_ (14 
SO = [ari ee 


where a=f/fo. This integral, when integrated between 
a, and de gives 


Si vonvar 


. Sd sa 1+a 
=B"fe| = ee log( "x 


l—a” l1-—a;, 1 —a, 





=] (15) 











Say ; a+1 
“( —a? * 8 EN] (16) 
= B” fol F(a1) — F(a2)]. (17) 


Table II gives values of a/1—a? and F(a) for several 
values of a. 











TABLE II. 
a a/(1—a*) F(a) a a/(1—a’) F(a) 
0.60 0.94 0.12 1.01 50.3 23.8 
70 1.37 .25 1.02 25.2 11.46 
75 1.71 37 1.03 16.9 7.41 
.80 2.22 56 1.04 12.7 5.39 
85 3.06 .90 1.05 10.2 4.19 
90 4.73 1.62 1.06 8.58 3.40 
.92 5.99 2.21 1.08 6.49 2.39 
.94 8.08 3.17 1.10 5.24 1.86 
95 9.75 4.00 1.15 3.57 1.12 
.96. 12.3 5.16 1.20 2.73 0.77 
.97 16.4 7.16 1.30 1.88 43 
.98 24.8 11.2 1.40 1.46 aa 
99 49.7 23.5 1.60 1.03 .14 
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The Effect of Nascent Hydrogen on the Magnetic Properties of Iron* 


R. K. REBER, Sloane Physics Laboratory, Yale University 
(Received July 9, 1934) 


Hydrogen was electrolyzed into flat ring cathodes of 
magnetically soft iron in 0.5 N solutions of H2SO, or of 
KOH. The design of the electrolytic cell permitted mag- 
netic testing at intervals. Magnetic hardening began after 
many hours provided the electrolytic current density was 
kept above a lower limit depending on the electrolyte. 
This hardening proceeded, at adequate current density, 
to an end point beyond which longer charging at any 
current density rarely produced any further changes in 


magnetic properties. Removal of hydrogen by storage in 
air at room temperature for 5 months, or by annealing at 
400°C in vacuum for 17 hours, did not diminish the 
magnetic hardness produced during electrolysis. It is 
suggested that the iron is locally cold-worked by non- 
uniform distribution in it of hydrogen in excess of the 
solubility limit at room temperature, and that uniformly 
distributed hydrogen up to this limit has very little effect 
on the ferromagnetic properties. 





HisTORICAL INTRODUCTION 


ASCENT hydrogen (H;*) liberated in 

contact with iron in an acid bath, or by 
electrolysis, rapidly embrittles ordinary fine- 
grained metal. Such brittleness is relatively 
transient, even at room temperature.' Magnetic 
effects of such treatments are not striking and 
have not been studied in detail.” 

The hardness of electrodeposited iron is 
usually ascribed to its retained hydrogen.’ 
Removal of hydrogen by heating is accompanied 
by magnetic, as well as mechanical, softening.* 
It seems possible, however, to get out the 
hydrogen without much altering the hardness,° 
and the structure of such metal is so anisotropic 
that interpretation of magnetic changes in it is 
complicated at best. We do not consider electro- 
deposited iron further herein. 

When this study was undertaken it seemed 


*Condensed from a Dissertation presented to the 
Faculty of the Graduate School of Yale University in 
Candidacy forsthe Degree of Doctor of Philosophy. 

1See e.g., L. B. Pfeil, Proc. Roy. Soc. A112, 182 (1926), 
and, most recently D. Alexejew, P. Afanassjew and W. 
Ostroumow, Zeits. f. Elektrochemie 40, 92 (1934). 

2D. E. Hughes, J. Soc. Telegraph Eng. 9, 170, 178 
(1880), detected no effect. H. Gries and H. Esser, Archiv f. 
Elektrot. 22, 145 (1929), found 3 mm rods of soft iron 
slightly hardened, magnetically, by 30 min. in 0.5 N H2SO, 
without current flow. J. Coulson, Phys. Rev. 19, 528; 20, 51 
(1922) observed that the magnetic moments of steel bar 
magnet cathodes in 25 percent H.SO, dropped by as much 
as 20 percent in some cases. He did not measure other 
magnetic properties. 

°L. Cailletet, Comptes Rendus 80, 319 (1875), showed 
hydrogen might be retained in large amount. 

*C. F. Burgess and C. Hambuechen, Trans. Am. 
Electrochem. Soc. 5, 201 (1904); R. Hugues, Rev. de 
métallurgie 22 [mém.], 764 (1925). 

® Guichard, Clausmann, Billon and Lanthony, Comptes 
rendus 192, 623, 1096; 193, 1084 (1931). 


possible that the extreme magnetic softness of 
iron treated in hydrogen at very high temper- 
atures by Cioffi® might, as he at first supposed, 
depend upon retained hydrogen. His later work’ 
showed that the process is one of purification. 
Sieverts® is probably right in supposing that 
excess hydrogen is rapidly given up by iron at 
any annealing temperature. Such considerations 
as these led us to decide that the effects of 
hydrogen on the ferromagnetism of iron would 
prove easiest to interpret if produced at room 
temperature in the purest metal we could prepare 
in pieces suitable for magnetic testing. This 
limited us to electrolytic means and suggested 
careful purification by Cioffiis method as a 
standard preliminary process. As it seemed likely 
that transient effects might occur, means had to 
be provided for making magnetic measurements 
during electrolysis. Mechanical strains had to be 
avoided, since these may entirely alter the mag- 
netic properties of soft magnetic materials. The 
saturation of iron with hydrogen at room tem- 
perature is a slow process.’ It seemed, therefore, 
that very long times of experiment might be 
necessary, and that pieces of different thickness 
had better be studied under otherwise identical 
conditions. This sharply limited the number of 
materials we could hope to investigate. 


*P. P. Cioffi, Nature 126, 200 (1930). 

7P. P. Cioffi, J. Frank. Inst. 212, 601 (1931); Phys. Rev. 
39, 363 (1932). 

3 A. Sieverts, Zeits. f. Metallkde. 21, 37 (1929). 

®T. S. Fuller, Trans. Am. Electrochem. Soc. 36, 113 
(1919), found, for example, that hydrogen liberated 
continuously on one side of a 1.6 mm iron wall began to 
escape from the other side after upwards of 15 hours, but 
that no steady state of flow had been reached in 200 hours. 
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298 R. K. REBER 
TABLE I. Experimental treatments. TABLE II. Magnetic properties—electrolysis. 
Current Rapid 
density change 
Thickness Solution amp. Ring Bmax B, H. between 
Ring Mate- em (10)73 (0.5 N) cm~? Time No. Before After Before After Before After hrs. 
No. rial Stock Ring of *(10)-3 hr. Notes 
1 23300 8550 13070 5300 0.227 0.332 3-15 
1 CRS 56 56 H2SO. 12 0-44 = (a) (b) 2 21400 6720 13800 5690 267.418 8-18 
24 44-73 (b) 3 16300 16300 9230 9230 182.182 _ 
48 73-121 (b) 4 14800 14800 13200 13150 .570 .574 — 
2 CRS 56 56 H2SOs 12 0-64 = (a) (c) 5 20400 11900 11190 6120 180.180 (?) 
3 CRS 56 56 H2SO. 3 0-240 (a) (d) 6 11900 2800 15700 6710 607 1.20 11-48 
4 CRS 56 56 KOH 160 0-120 (a) 7 32400 6240 12140 3840 135 = .230 0-55 
400 120-216 8 21900 7750 9800 4400 -151 .257 60-85 
5 CRS 56 56 KOH 1400 0-48 = (a) (d) 9 19000 20300 13500 12100 243.225 — 
6 CRS 56 56 H2SO4 12 0-96 (e) 10 35800 15860 10440 4420 211 -246 16-110 
7 CRS 216 216 H2SO. 12 0-104 (f) (g) 11 29900 13800 12660 4930 157.168 92-210 
8 CRS 216 216 KOH 12 0-48 (f) 12 25130 12430 10670 5060 -148 .210 96-154 
H2SO.4 12 48-169 (c) (h) 12430 4720 5060 3070 210 =.391 203-300 
9 CRS 216 216 H2SO,4 0 0-200 (f) (i) 13 19300 4600 7410 3320 -138 -380 95-123 
10 CRS 317 317 H2SO« 12 0-115 ) (j) 14 18000 6080 6280 2610 170.255 66-96 
11 CRS 317. 317 H2SO.4 12 0-228 15 17100 17000 5820 5790 -151 -150 — 
12 Armco 79 79 H2SO,4 12 0-203 te) (k) 16 15960 2460 4120 1980 -103 475 155-224 
48 203-231 17 15800 15800 8000 8000 184.184 — 
3 231-252 18 20200 3970 9320 3000 142.449 165-195 
48 252-324 
13 Armco 79 79 H2SO,4 48 0-161 (a) 
14 Armco 79 79 H2SO4 192 0-60 (a) 
0 60-62 ° ° . . 
192 62-120 not attain maximum possible purity, partly, no 
15 Armco 79 79 H2SO« 6 0-336 (a) (d) . 7 
16 Armco 316 316 H:SOs 8 0-330) () doubt, because the hydrogen supplied to our 
192 330-374 . 
17 CRS 316 56  HSOs 12 0-100 (f)(m)(n) furnace passed over hot alundum before reaching 
18 Armco 316 56 H2SO,4 48 0-208 (f) (m) A . . 
the specimens. All rings were magnetically soft 
(a) Heat treatment: 6 hours at 1500°C in H2 at atmospheric pressure; enough, nevertheless, to make even slight hard- 


cooled in He to 880°C; pumped down to a pressure of the order 
of 0.1 mm; 2 hoursat 880°C; furnace cooled to room temperature. 

(b) Electrolysis stopped during magne tic tests. This allowed too much 
attack by acid so that in all other cases (except in Ring No. 12) 
the elec troly sis was continued at very low current density during 
magnetic tests. 

(c) Magnetic field in ring (H =18) throughout, except as necessarily 
removed for magnetic testing. 

(d) Electrolysis in an open tank (no magnetic tests during progress). 

(e) High temperature treatment in hydrogen omitted. Annealed 3 
hours at 850°C in vacuum. 

(f) Heat tre atment: same as (a) except for longer times, 
1500°C and 3 hours at 880°C. 

(g) Subsequently placed in a toroidal wooden box and wound with 
primary and secondary. Magnetically tested at intervals for 5 
months. See Table III. 

(h) Subsequently annealed at various temperatures. See Table IV. 

(i) Ring so much attac ked by acid that magnetic measurements were 


12 hours at 


no longer possible even with corrections for change in cross 
section. 

(j) Charging current interrupted for an unknown period not exceeding 
8 hours. 

(k) In dilute HNOs one hour before placing in electrolytic cell. Elec- 
trolytic current not turned down during magnetic testing. 


Results of magnetic tests corrected for presence of electrolytic 
current in a separate experiment. 

(1) In dilute H2SO« one hour before placing in electrolytic cell. 

(m) Ground down before anneal at 880°C. 

(n) Piece fell out of ring, preventing furthe r study of magnetic proper- 
ties. 


EXPERIMENTAL METHODS AND RESULTS 


The specimens were rings (6.35 cm outside 
diameter, 5.08 cm inside diameter) punched or 
lathe cut from sheets of cold rolled steel (CRS 
in Table I) or of Armco iron. All rings of a given 
stock thickness came from a single sheet. The 
standard treatments before use, Table I, notes 
(a), (f), should have reduced ordinary impurities 
to the point where chemical analysis is unin- 
formative, and analysis was therefore not at- 
tempted.” Magnetic results show that we did 


1° We are indebted to Mr. Cioffi for information regarding 
his methods and results prior to publication. 


ening easily apparent. We attach no importance 
to the absolute values of the magnetic constants 
as reported in Table II and, except as noted in 
the discussion, data for one ring are to be com- 
pared only with earlier or later data for the same 
ring. 

Fig. 1 shows the electrolytic cell i in which rings 
were charged with hydrogen under conditions 
described in Table I, and in which magnetic 
characteristics presented in Table II were ob- 
tained. The solution was kept flowing through 
the cell throughout a test, and loss of water was 
made good at frequent intervals. Charging 
current density was generally reduced below 
(10)-* amp. cm during magnetic tests to 
prevent “cross-talk” between charging and gal- 
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Fic. 1. Electrolytic cell (sections). 





























MAGNETIC PROPERTIES OF IRON 299 
TABLE III. Magnetic properties—aging at room temperature. TABLE IV. Magnetic properties—annealing. 
Ring No. 7 jimax B,; H, Ring No. 8 lan B, H, 
Before electrolysis 32400 12140 0.135 Before electrolysis 21900 9800 0.151 
0 hr. after electrolysis 6250 3710 .236 After electrolysis 7750 4400 .257 
14 hr. after electrolysis 6250 3770 .235 = After 17 hr. at 400°C 7800 4400 .256 
40 hr. after electrolysis 6250 3790 .233. After 12 hr. at 550°C 9750 6350 aan 
13 days after electrolysis 6250 3810 .233 After 12 hr. at 830°C 15800 7650 .180 
3 months after electrolysis 6240 3840 .231 
5 months after electrolysis 6240 3840 .230 








vanometer circuits. The rings fitted the cell so 
loosely as to prevent strain in assembly or as a 
result of ordinary temperature changes. Charging 
from both sides (in a vertical plane) tended to 
reduce mechanical and magnetic asymmetry. 
No change in dimensions of as much as one part 
in 10,000 could be detected in experiments (not 
reported in Table I) where such changes were 
looked for. 

Magnetic measurements were made by bal- 
listic galvanometer (Leeds and Northrup Type 
R, sensitivity 0.0003 microcoulomb per scale 
division). A typical winding (not shown in Fig. 1) 
might consist of a primary of 300 turns of No. 26 
(0.0405 cm) enameled copper wire, with a 
secondary of 1000 turns of No. 32 (0.0202 cm) 
enamelled copper wire, both distributed over the 
whole circumference. Complete magnetization 
curves and hysteresis loops, usually up to 
Hyax = 45, were taken only at beginning and end 
of an electrolysis. During a run only enough 
points were taken to give approximate values of 
max, B, and H, at chosen epochs. An effort was 
made to fix the times at which any important 
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Fic. 2. Magnetization curves for ring No. 16. 


change began and ended. Table II is a condensed 
account of these major effects. 

Figs. 2, 3 and 4 show typical magnetic results. 
The most striking feature is that even when 
magnetic properties eventually change most 
profoundly there may be little or no change for 
a long time, then a relatively rapid change, and 
finally another period of constancy. All large 
and delayed changes of this sort increase mag- 
netic hardness. If there is any immediate effect 
it is a relatively small decrease in magnetic 
hardness (Fig. 4). 

The results of a few special aging and an- 
nealing tests are presented in Tables III and IV. 
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Fic. 3. Hysteresis loops for ring No. 12. 
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Fic. 4. umax, Br, H- as functions of time of electrolysis for 
ring No. 2. 
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DISCUSSION 


These experiments and others available for 
comparison seem consistent with the following 
simple picture of what happens when hydrogen 
is liberated at the surface of an iron specimen in 
an electrolytic cell. 

Mere immersion in dilute acid, without any 
current flowing, relieves surface strains and 
produces slight magnetic softening. This change 
is not observed in KOH solutions, is relatively 
more important in thinner rings, and can be 
prevented from occurring during electrolysis in 
H.SO, solutions by preliminary etching with 
HNOs, which does not liberate H*+. Small ad- 
herent grains of alundum fall off during this 
initial softening or during any treatment which 
prevents its occurrence. 

Nascent hydrogen enters iron readily at grain 
boundaries or other crystal imperfections which 
reach the surface of the piece. Penetration of 
relatively thick pieces takes but a short time— 
of the order of one hour or less—results in me- 
chanical brittleness in fine-grained metal," but 
has almost no magnetic effects. In the present 
investigation the grains had dimensions of the 
order of a millimeter or more and embrittlement 
was not observed. 

Hydrogen diffuses more slowly into the more 
nearly perfect parts of the grains and there is 
some evidence for the idea that the larger such 
regions are the longer it takes to saturate the 
entire specimen with hydrogen, supplied at the 
same adequate rate in all cases. At least it was 
observed here that metal of greater purity, as 
judged by larger grain size and higher corrosion 
resistance, required a longer time for the onset 
of profound changes of magnetic behavior. Dif- 
fusion of hydrogen in iron is probably atomic 
rather than molecular, to judge by its dependence 
on current density.” The relatively small effect 
of thickness upon the time required for saturation 
(onset of magnetic changes) suggests that the 
diffusion coefficient increases as saturation is 
approached. Rings 16 and 18 may be cited in 
support of this idea. The fact that abrupt mag- 








1 J. H. Andrew, Trans. Faraday Soc. 9, 316 (1913); 14, 
232 (1919). 

12 Probably proportional to the square root of current 
density. See M. Bodenstein, Zeits. f. Elektrochemie 28, 517 
(1922); G. Borelius and S. Lindblom, Ann. d. Physik 82, 
201 (1927). 
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netic changes still occur with relatively high 
current density also favors this hypothesis. 

If we regard the iron as just saturated with 
hydrogen when important magnetic changes 
begin it follows at once that isotropic annealed 
iron saturated with hydrogen at room temper- 
ature is still magnetically soft. Magnetic harden- 
ing must then be associated with some sort of 
“supersaturation.’’ We propose to use this term 
with some caution, for the easy escape of hy- 
drogen from embrittled iron shows that at least 
a part of the metal may already be supersatu- 
rated at a very early stage of the process. It has 
been suggested by Alexejew and co-workers! that 
the hydrogen thus loosely bound is concentrated 
in grain boundaries or in ultramicroscopic cracks. 
We must conclude, if this is correct, that hy- 
drogen may escape from the free surface of an 
iron wall befovre all the grains composing the wall 
have been supersaturated in the sense here 
intended. 

There seems to be a threshold current density, 
of the order of (10)-? amp. cm~ in the 0.5N 
H2SO, here used, below which large magnetic 
effects are not observed. In KOH solutions at the 
same normality this threshold is much higher, 
and it probably depends in both cases upon the 
purity of the electrolyte, for Alexejew! found no 
embrittlement in HeSO, solutions free from H2S 
and other hydride forming impurities. Whether 
this threshold is exceeded little or much, about 
the same magnetic hardening finally occurs, the 
only important differences lying in the time for 
onset of the change and in the period for its 
completion. These times are of the order of those 
apparently necessary for the attainment of 
stationary conditions in the transport of hy- 
drogen through iron sheets of comparable thick- 
nesses, as in the experiments of Fuller? and 
others.'* 

The magnetic results, when notable change 
occurs, suggest that the iron has been cold 
worked. Not only do comparable changes result 
from over-strain by stretching or rolling, but the 
low temperature for recovery—below 550°C—is 
characteristic of plastically deformed iron. We 
find a plausible process for such cold work in the 
internal liberation of hydrogen from super- 


13 Most recently commented upon by C. A. Edwards, 
Nature 133, 379 (1934). 











saturated metal. The threshold current density, 
on this view, is that necessary to maintain flow 
through the iron against a back pressure great 
enough to produce plastic deformation. At a 
given current density the volume of cold-worked 
metal rises rather rapidly as soon as the process 
begins and the process stops as soon as all the 
metal is strain-hardened. Tests show that the 
entire volume of the rings is affected and that 
the condition is unaltered by low temperature 
processes which remove hydrogen but do not 
promote recrystallization. 

During the cold-working process, at least, the 
hydrogen is non-uniformly distributed, which 
explains why so low an average concentration, 
probably under one atom of hydrogen per 
hundred atoms of iron,“ can be effective. The 
small difference in the time for complete harden- 
ing in rings of different thicknesses makes it 
probable that hydrogen is being fed at once to 
all points of weakness through many channels of 
low resistance rather than entering by plane 
diffusion from the exposed surface. Coulson’s 
experiments? on permanent magnet steels show 
particularly rapid changes. 

Cooling iron saturated at a high temperature— 
in the annealing range or higher—may result in 
temporary brittleness but no permanent hard- 


4 R. Chandler, J. Soc. Telegraph Eng. 9, 168 (1880); M. 
Thoma, Zeits. f. physik. Chemie 3, 69 (1889); A. Thiel and 
W. Hammerschmidt, Zeits. f. anorg. Chemie 132, 15 
(1923); H. QO. von Samson-Himmelstjerna, Zeits. f. 
anorg. Chemie 186, 337 (1930). 

%E. Heyn, Stahl u. Eisen 21, 913 (1901); Metallo- 
graphist 6, 39 (1903); H. Jellinek, Zeits. f. Physik 66, 543 
(1930). 
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ening has been detected by either mechanical or 
magnetic tests. 

We must suppose that embrittlement by the 
initial entry of hydrogen is not accompanied by 
any cold work of the sort here postulated, for if 
it were it is hard to see why it should not im- 
mediately progress at the expense of soft mag- 
netic material and give a much earlier onset of 
magnetic hardening. The sharp separation be- 
tween magnetically ineffective saturation and 
magnetically effective supersaturation must be 
insisted upon. 

Bridgman’® has reported premature and pro- 
gressive deterioration of steel under very high 
pressures of molecular hydrogen at room tem- 
perature. This suggests that the cold work we 
postulate may occur under any conditions which 
supersaturate part of the metal, especially if 
high elastic strains already exist. 


CONCLUSION 


Nothing has been said about any effects of 
magnetic fields superposed during electrolysis 
(rings 2 and 8), nor about the possible accelera- 
tion of magnetic changes by corrosion. These and 
other points of interest must be left unsettled. 
It is hoped that this publication will call atten- 
tion to the imperfections of the theory outlined 
above, and stimulate its correction or replace- 
ment. 

The author is indebted to Professor L. W. 


-McKeehan for suggesting this investigation and 
for help in arranging the results for publication. 


16 P, W. Bridgman, Proc. Am. Acad. 59, 173 (1924). 
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The Piezoelectric Properties of Quartz and Tourmaline 


GERALD W. Fox AND GEorGE A. Fink, Physics Laboratory, Iowa State College 
(Received January 23, 1934) 


The piezoelectric properties of quartz and tourmaline have been investigated through a 
study of the converse piezoelectric effect. In the case of quartz, the piezoelectric constant has 
been measured from point to point across the specimen. Considerable variation in the ‘‘con- 
stant’’ has been observed in different samples of quartz. Attention is called to difference 
between the motions of piezoelectric plates when subjected to static electric fields and when 


subjected to alternating electric fields. 





INTRODUCTION, 


N preparing quartz piezoelectric plates, it 
frequently occurs that blanks cut from neigh- 

boring positions in a crystal of optical quality, 
even though treated identically in the finishing 
process, produce oscillators of widely different 
abilities. Methods are available for locally testing 
oscillator blanks by applying pressure and ob- 
serving the resultant charge developed (direct 
piezoelectric effect) and, although these show 
that there exists considerable non-uniformity in 
quartz, they give little information about the 
absolute value of the piezoelectric constant from 
point to point. It was, therefore, thought worth 
while to investigate the local variations in the 
piezoelectric quality of quartz by using the 
converse effect wherein an electric field is applied 
to the plate and the change in thickness meas- 
ured. 

Most previous determinations of piezoelectric 
constants by the direct effect have necessarily 
given only the average for the specimen tested 
since no means was provided for determining 
how much of the charge liberated came from 
any particular area of the face.':?*3 Even a 
method using concentrated loads, such as that 
of Dawson,‘ is questionable for evaluating a 
piezoelectric constant at a point. A contact 
point with dimensions smaller than the thickness 
of the plate being tested sets up a very compli- 
cated system of stresses and strains. Under 
these conditions, it seems rather doubtful 
whether the charge collected by the contact 


! J. Curie and P. Curie, J. de physique [2] 1, 246 (1882). 

? E. Riecke, and W. Voigt, Ann. d. Physik 45, 549 (1892). 

?W. C. Réntgen, and A. Joffe, Ann. d. Physik [4] 41, 
490 (1913). 

* Dawson, Phys. Rev. 29, 553 (1927). 


point or probe is simply equal to the product of 
the applied force and the piezoelectric constant. 
In the first place, the contact point or probe 
may not collect all the charge produced if the 
quartz does not have a conducting coat, since 
portions of the quartz surface, not in contact 
with the probe, are strained and produce surface 
charges. Secondly, the proportionality between 
strain and polarization may fail when the stress 
producing them is very high. Thirdly, at many 
points in the quartz, the stress is not acting 
along the normal to the surfaces and other 
constants play a part beside the one appropriate 
to uniform stress along the normal to the surfaces 
of the quartz plate. 

Even in the present method, elements of the 
plate are not free to expand independently but 
are elastically bound together. However, the 
applied field is uniform. To get a corresponding 
state of affairs with the direct effect, it would be 
necessary to apply uniform pressure over the 
faces of the quartz and measure the charges 
produced on selected areas of the surface. The 
only previous determination by the converse 


_effect known to the authors, that of (Ny Tsi)Ze,' 


measured an average value for the specimen. 
Ze’s method gave the interesting result of an 
apparent saturation with strong fields, but he 
could only test a long thin piece or a stack of 
plates. 

To test a single piece the size of an oscillator 
blank by the converse effect, an extremely 
sensitive method is necessary for measuring the 
resulting thickness changes. These are only of 
the order of 10-7 cm at an applied potential 


5 Ny Tsi Ze, J. de physique et le radium [6] 9, Part 1, 
page 13 (1928). 
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difference of 500 volts. The amplifying lever 
developed for these measurements has already 
been described.® 


RESULTS 


Variation of constant from point to point 


This lever and connected apparatus, which 
might be called a piezoelectric dilatometer, was 
first set up with three foil-coated plates of glass 
about two millimeters thick under its three knife 
edges. No movement of the interference fringes 
used for observation could be seen when potential 
differences up to 500 volts were applied to the 
plates. A movement of one-tenth of a fringe, 
corresponding to a change in thickness of 0.15A, 
could have been detected. The fringe shifts 
observed later with average quartz plates were 
in the neighborhood of six or seven times the 
distance between fringes. A quartz plate with 
faces normal to the optic axis and a piece of 
fused quartz gave shifts of 1.2 and 1.3 fringes 
on 500 volts when substituted for the middle 
glass plate, but this, being the same for both 
directions of the applied field, was an electro- 
static or electrostrictive effect. 

Three blanks, numbers 37, 38 and 39, that 
had been cut with faces normal to the three 
electric axes of the same crystal, were tested and 
the positive directions of the axes were found to 
agree with Fig. 1. The positive direction of an 
electric axis is the direction of the electric field 
which produces an expansion of the crystal 
along the same axis. The piezoelectric deforma- 
tions were nearly proportional to the applied 
voltages, as may be seen from the data plotted 
in Fig. 2. The points are not quite on a straight 
line; showing that, besides the linear piezoelectric 
effect, there is an electrostrictive or other effect 
proportional to the square of the voltage. If the 
average of two fringe shifts is taken, one with 
the field equal in strength but opposite in 
direction, those effects proportional to the square 
of the voltage will be eliminated. This has been 
done to obtain the points in Fig. 3 which are 
seen to lie on a straight line to within the limit 
of experiment error. The crosses in Figs. 2 and 3 
represent readings taken with the center knife 
edge of the dilatometer resting on a spot near 





® Fox and Fink, Rev. Sci. Inst. 4, 276 (1933). 
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Fic. 1. End view of quartz crystal showing X- and Y-cuts. 
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Fic, 2. Data for plate 37 with varying voltage. 
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Fic. 3. Data of Fig. 2 averaged for + and — voltages. 


the center of plate 37. The plate dimensions 
were about 12152 mm. The piezoelectric 
constant, dy in the usual notation,’ calculated 
from these data, is 5.17A/e.s.u. With the plate 


7W. G. Cady, Int. Crit. Tab. 6, 208. 
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moved a few millimeters so that the knife edge 
rested on a different part of the plate, readings 
were obtained giving dy, equal to 6.04A/e.s.u. 
These readings are represented by the dots in 
Figs. 2 and 3. 

The question then arose whether the difference 
is due to an actual difference in the piezoelectric 
properties of the quartz at the two points or to 
changes in the apparatus caused by moving the 
quartz and readjusting the amplifying lever. 
To answer this question, another test was made 
at the point where the value of 6.04A/e.s.u. was 
obtained. With the dilatometer readjusted after 
standing idle a day or two, the resvlt was 
6.20A/e.s.u. The readings for this run are repre- 
sented by the circles in Figs. 2 and 3. Next, the 
specimen was turned over and readjusted again, 
taking care to have the knife edge resting as 
nearly as possible on the same part of the 
specimen. 5.97A/e.s.u. was the constant found 
this time. The data are plotted with triangles in 
Figs. 2 and 3. The variation in the three values 
6.04, 6.20 and 5.97 is about the limit of experi- 
mental error. However, the difference between 
the average value at the second spot and the 
value obtained at the first spot is definitely 
beyond the limit of experimental error. Two 
other X-cut plates from the same crystal gave 
the values 6.1 and 5.5A/e.s.u. at one spot on each. 

Three Y-cut blanks from the same crystal 
were tested and the changes in thickness were 
found to be much smaller than for the X-cut 
blanks. Values for the piezoelectric constant at 
one point on each were: 1.4, 1.1 and 2.2A/e.s.u. 
The constant connecting the applied voltage 
and the resulting change in thickness in the 
case of a Y-cut plate is not dy but it is not clear 
just what it should be. According to the ordinary 
theory, assuming uniform structure and equal 
“strengths” of the three axes at every point, 
the Y-cut plates should undergo pure shear with 
no change in thickness. The changes in thickness 
that were found for this type of plate indicate 
that the three electric axes at a point in a piece 
of quartz do not have the same strengths. The 
orientation of the plates with respect to the axes 
of the crystal was accurate to within a degree. 
Errors in orientation would then be expected to 
produce an expansion of the Y-cut plates less 
than three percent of the expansion of the X-cut 
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Fic. 4. Deformation of Y-cut plate in an electric field. 


plates; i.e., less than 0.3A/e.s.u., which is much 
less than the values found. 

A possible explanation on the basis of unequal 
strengths is illustrated in Fig. 4. The positive 
directions of the electric axes in a piece of quartz 
are shown by therarrows 1, 2 and 3. If the field 
is in the direction shown, it will have no compo- 
nent in the direction of axis 1. It will have a 
positive component in the direction of axis 2 
and will, therefore, cause expansion in that 
direction. In the direction of axis 3, the field 
component is negative and, therefore, will 
produce contraction. If the elongation per unit 
length in the direction of axis 2 or AD is equal to 
the contraction per unit length in the direction of 
BC, then the rectangle A BDC will undergo pure 
shear and become the parallelogram ABD’C’ 
with no change in altitude, at least to a first 
approximation. However, if the elongation and 
contraction per unit length are not equal, the 
deformation will have some other nature than 
that of a pure shear and the resulting parallelo- 
gram will have an altitude differing from that 
of the rectangle. It is also possible that some of 
the apparent changes in thickness of the Y-cut 
plates—and the X-cut plates, too—may be due 
to bending of the plates. Bending would be 
caused by local variations in the expansion along 
an electric axis and by variations in the corre- 
sponding lateral contraction at right angles to 
this axis. 

A specimen of tourmaline 2.39 mm thick was 
tested next. It was cut with its faces normal to 
the principal axis which is both the optic and 
electric axis in this crystal. The piece from 
which the specimen was cut was roughly tri- 
angular in cross section, about 8 mm on a side, 
green in color except in the center which was 
yellow. The value of 5.2A/e.s.u. was found at 
one spot and 5.1A/e.s.u. at approximately the 
same place with the other side up. The values 
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Fic. 5. Variation of the piezoelectric constant across a plate. 


of d3; for tourmaline given by Cady’ are 5.78, 
5.4, and 5.6A/e.s.u. 

After these first tests indicated that the piezo- 
electric “‘constant” varied from point to point 
across the face of a specimen, the apparatus was 
changed so that specimens could be moved under 
the center knife edge without removing the 
protecting covering and taking the apparatus 
apart. With this means for moving the specimen, 
quartz plate number 37 was retested at twelve 
points approximately a millimeter apart in a 
line across its face. The results are shown by 
circles in Fig. 5. Upon going across the face 
again, testing approximately the same points, 
the results were quite different as shown by the 
crosses in Fig. 5. The knife edge which rested 
on a narrow strip of the specimen’s foil coating 
less than 0.01 mm wide would probably not 
touch exactly the same spots the second time, 
but the variations are probably less than } mm, 
certainly much less than 2 mm, the thickness 
of the plate. It seems unlikely that the change 
in thickness of a plate 2 mm thick would vary 
much from one point on its surface to another a 
half mm away or less, since, to do so, would set 
up concentrated shearing and bending stresses. 


Variation of constants with surface conditions 
and repeated application of the electric field 


The wide discrepancy between the values 
shown in Fig. 5 for two traverses across the same 
quartz plate may be due to some of the following 
causes ; different surface conditions, the effects of 
repeated application of the electric field, possibly 
changes in the contact between the plate and 
its support as the plate slides over its support. 
In any event, further investigations of these 
peculiar results must be made before an adequate 
explanation can be offered. 
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Field Reversed 


Fic. 6. Showing reversal of fringe movement with reversal 
of field. 


In some further experiments with quartz and 
glass plates, effects due to repeated application 
of the electric field were observed. This was 
especially true of plates that had been coated 
by cathode sputtering instead of by the original 
method of coating in which aluminum foil was 
stuck on with a very small quantity of oil. The 
movement of the fringes observed upon the first 
application of the electric field was usually not 
as great as that attained by repeated application. 
In one experiment with a glass plate having a 
sputtered coat, the response when the field was 
applied the fifteenth time at intervals of about 
a half minute was nearly three times as great as 
that obtained the first time. Still more interesting 
was the behavior when the field was reversed 
for the first time. The first movement of the 
fringes, when the field was applied in the reversed 
direction, was in the same direction as that 
produced by the original direction of field. This 
movement would slow down, stop, and reverse, 
carrying the fringes back to their zero position 
and beyond to a final steady deflection opposite 
to that with the original field. The movement 
out and back to zero took ten or fifteen seconds, 
while a steady state would be reached in twenty 
to thirty seconds. When the field was applied 
again in the second direction, the reversal of 
movement did not occur, the fringes moving 
immediately in the direction of the final deflec- 
tion. Some typical results are shown in Fig. 6. 
The final ‘“‘saturated”’ fringe movements are seen 
to be fairly consistent. 

At some times, a time lag of the fringe move- 
ment was observed greater than that due to 
damping of the motion of the amplifying lever. 
Records of this lag were made by recording on 
a chronograph the times of passage of successive 
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Fic. 7. Curves showing time lag of fringe movements. 


fringes across the intersection of the cross-hairs 
in the observing telescope. Fig. 7 contains data 
from two such runs. The slight slope of the flat 
tops of the curves may be due to leakage currents 
heating the specimen. 


Comparison of an excellent with a worthless 
blank 

The most satisfying results were obtained from 
a comparison of a specimen from a_ badly 
twinned and flawed quartz crystal with a plate 
that oscillated even though its edges were not 
carefully ground as is necessary with some 
inferior blanks. The data plotted in Fig. 8 show 
that the quartz of poor quality showed a much 
smaller response at every point than the good 
quartz, even reversing its sign at one point. 
The good oscillator gave a large response at 
every point tested and, except for the lowest 
point, all are higher than any obtained with any 
other specimen and higher than the values 
reported by previous observers using the direct 
effect. 


REMARKS ON MOTIONS OF PLATES UNDER 
RESONANCE AND STATIC CONDITIONS 


It is well to mention that the motions of 
quartz plates observed in this investigation are 
not necessarily the ones that the same plates 
would make in the oscillating circuits where 
they are customarily used. Though there is no 
doubt that when a static electric field is applied 
to an X-cut plate, as was done in this investi- 
gation, there is a decided expansion or con- 
traction in the direction of the field, there is 
considerable doubt that the major surfaces of 
the plate vibrate as a whole, as is customarily 
supposed, when the plate is driven at some 
resonance frequency. The experiments of Oster- 


> G&G. A. FIRES 
lo> 4° 
4 Agsu . 
8-4 e e . ° 
6- . 
4- Good 
| Non- Oscillator Oscillator 
2] #*4) . #4 
fs) , : _ 





-2 TYrTT rr TTT T 
o ¢ wb. # 
mm from end of plate 


Fic. 8. Comparison of piezoelectric quality of good and 
poor quartz. 








ee ait TTTTT 
/ 4s 


berg®: * 1% "! show that the actual motions of 
quartz plates driven at resonance frequencies 
are very complex, and are frequently due to 
coupling between waves reflected from the faces 
and edges of the plate. The actual motion of the 
surfaces of a plate is then the resultant of a 
number of periodic disturbances propagated 
through the plate. Koga™ has recently developed 
a theory which states that the longitudinal 
thickness vibration is to be expected for the 
X-cut plates. This has usually been assumed to 
be the mode of vibration for this type of plate. 
Interferometer studies, however, have not borne 
out his idea. The motions of the surfaces of 
such a plate seem closely dependent upon the 
shape of the plate and how accurately it has 
been cut with respect to the faces of the crystal 
but in no case has actual motion of the major 
surfaces as a whole been observed except in 


recent studies by Straubel who states that, by 


a suitable choice of plate shape, it is possible to 
produce oscillators having a single resonance 
frequency in which whole faces vibrate as a 
unit according to his powder figure studies. 
According to Koga’s theory, the so-called parallel 
cut always vibrates in the pure shear mode but, 
here again, interferometer studies do not bear 
this out."* In the static case for the Y-cut, 
shearing does exist just as expansion along the 
electric axis exists in the X-cut but it should 
not be assumed that the motions of a plate under 
static electric fields are the same as those it would 


execute if vibrating at a resonance frequency. 


8H. Osterberg, Proc. Nat. Acad. Sci. 15, 892 (1929). 
°H. Osterberg, J. Opt. Soc. Am. 22, 19 (1932). 

10H. Osterberg, J. Opt. Soc. Am. 23, 30 (1932). 
11H. Osterberg, Phys. Rev. 43, 819 (1933). 

121. Koga, Phil. Mag. (London) 16, 275 (1933). 

18 Straubel, Phys. Zeits. 32, 586; 32, 222 (1931). 

14 Private communications with Dr. Osterberg. 
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It is shown that the variation with depth of the density and elastic constants of a half-space 
can be determined from the observed motion of the surface caused by a given distribution of 


stresses on the surface. 





INTRODUCTION 


HE information one obtains from the inter- 

pretation, under the ray concept, of travel- 
time curves due to a seismic pulse is the variation 
with depth of the quantities n/p and (A+2y)/p. 
This method makes use only of the arrival-times 
‘of the P- and S-waves, but not of their form. 
The question then arises as to what additional 
information, if any, can be obtained from an 
analysis of the forms of these waves. Clearly, 
the relevant observational data are the actual 
motions of the surface of the earth during the 
entire passage of the shock, and not the seismo- 
gram itself. Let us assume then that from a 
knowledge of the constants of the seismographs 
we have reproduced in each case the actual 
motion of the earth. The resulting record at a 
given station will depend first of all on the time- 
and space-variation of the initial pulse. Assuming 
a definite form for the pulse, we can investigate 
the resulting motion of the surface when the 
physical properties of the medium are known. 
Conversely, if the form of the pulse and the 
motion of the surface are given, we may attempt 
to find the physical constants of the medium. 
We should expect that such an analysis would 
yield at least the information obtainable from 
the arrival-times alone. 

The problem at hand may be called, in 
Slichter’s terminology, an inverse boundary value 
problem in elasticity. Similar inverse boundary 
value problems in electrostatics and electro- 
dynamics have recently been treated by Slichter,! 
Langer? and Stevenson.’ In the following we shall 
discuss the general features of the forced vibra- 
tions of an heterogeneous elastic half-space and 


* Read at the 1934 Spring meeting of the Eastern 
Section of the American Seismological Society. 

1L. B. Slichter, Physics 3, 273 (1932); 4, 411 (1933). 

?R. E. Langer, Bull. Am. ‘Math. Soc. 39, 184 (1933). 

3A. F. Stevenson, Physics 5, 114 (1934). 


shall give a formal solution of the inverse 
problem for two relevant types of stresses 
applied at the surface. In the first case the 
dilatation is zero throughout and the analysis 
yields the density p and the elastic constant yp as 
functions of depth. In the second case the 
dilatation is finite and we obtain p, uw and the 
elastic constant \ as functions of depth. 


1. STATEMENT OF THE PROBLEM 


On the surface of a half-space whose elastic 
constants } and uw and density p vary only with 
depth is applied a given distribution of stresses; 
the functions d, u and p are to be determined from 
the observed motion of the surface. Let the axis of 
z point downward; then X, uw and p are functions 
of z only. The general equations of elastic motion 
in an heterogeneous solid : 
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reduce under these assumptions to 
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O7u 06 du fou dw 
p—= (A+p) —+aY?u+— (—+—), (1.5) 
of Ox dz \dz 0x 


0*v 06 du fdv dw 
p—=(A+nu) —+ wv? +— (tn) (1.6) 
of oy z\dz dy 


~ 


Pati 06 dd dypdw . 
= (A+) —— w+s—+2— —, (1.7) 
at dz dz dz 
Now assume that 
u=0F/dx+0G/dy, v=0F/dy—0dG/dx, 
w=dF/dz+w’, (1.8) 


where F, G, and w’ are functions of x, y, 2 and ¢. 
We have merely effected here a linear transfor- 
mation of the variables u, v and w into the 
variables F, G and w’. On substituting (1.8) in 
(1.5), (1.6) and (1.7) it is found that the resulting 
equations are satisfied by the following system 
of equations: 
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The boundary conditions are that the functions 
shall be bounded (except, perhaps, at points of 
discontinuity of the applied stresses) and that 
the stresses 


Pez =b+ 2y0w/d2, (1.12) 
= p(dw/dx+du/dz), (1.13) 
= u(dw/dy+dv/dz) (1.14) 


at the surface (z=0) shall equal the prescribed 
applied stresses. Noticing that 


5=V?F+0w’ /dz (1.15) 


PEKERIS 





it is clear from the above equations that also in 
a heterogeneous medium two types of elastic 
motion are possible:* one without dilatation 
corresponding to 


G+0, F=0, w’' =0, (1.16) 
and one with dilatation corresponding to 
G=0, F #0, w' £0. (1.17) 


Moreover, in the first case the motion is purely 
horizontal and the normal stress p., is every- 
where zero, while in the second case the vertical 
component of the curl of the displacement is 
everywhere zero. 


2. CASE OF A Non-DILATATIONAL DISTURBANCE 


Let us begin with a discussion of the solutions 
of (1.9) which are periodic in time: 
G(x, y, 2, t) =e'?*H(x, y, 2). (2.1) 


In a cylindrical system of coordinates, equation 
(1.9) becomes 
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The radial and azymuthal displacements, which 
we shall denote by u and v, respectively, are 
according to (1.8) given by 


u=(1/r)dG/08, v=—d0G/dr; (2.3) 
the corresponding stresses are seen to be 
=p(du/dz), (2.4) p:o=pu(dv/dz). (2.5) 


We shall limit ourselves to solutions which are 
symmetrical about the z-axis. Thus: 


u=p.,=0. (2.6) 
Under this assumption Eq. (2.2) is satisfied by 
H(r, 2) =R(r)-Z(z) 


@R 1dR 
(2.7) 


and 


4 This"result was given first by J. H. Jeans, Proc. Roy. 
Soc. A102, 554 (1923). 
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dz? 
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1 du dZ 
+- — ——(k?— p*p/u)Z =0. (2.8) 
uw dz dz 
An appropriate solution of (2.7) is the Bessel 
function Jo(kr). Of the solutions of (2.8), we 
adopt the one which for a fixed value of p and 





o(s, 7, t)= f eirtd p f o(p, RI (kr) Zi (2, k, p)dk. 


The function ¢(, k) has to be chosen so that 


(Pa) o-0= Hof endo | o(p, k)kIi(kr)Zi(0, k, p)dk=f(r, t), 
—oo 0 
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z(+0) is asymptotically like e~** for large values 
of k, and call it Z,(z). It follows now from (2.3) 
that 


f(r, t) being the stress applied at the surface. The quantity that is observed is 


v(0, nd=f cindy [ o(p, k)kJi(kr)Z,(0, k, p)dk. 
—o 0 


v=kJ,(kr)Z,(z)e'. (2.9) 
Let us seek solutions for v of the form 
(2.10) 
(2.11) 
(2.12) 


Now the integrals in (2.11) and (2.12) are of the type which can be inverted by an application of 


the Fourier identity: 


F()= (1/2) eindp | F(d)e-?\dn, 


and the Fourier-Bessel identity : 


F(r)= f “T.(kdkdk f ” F(R) Jo(RR)RAR, 
0 0 


(2.13) 


(2.14) 


© ioe) 
provided that f | F(t) | dt and f F(r)(r)‘dr exist and are absolutely convergent. We can therefore 
0 


—o 


write formally 


o(p, k)Z,(0, k, p)= (1/2n) f Terre [ v(0, r, t)e~*?dt, 
0 


(Pp, k)Z,(0, k, p)= (1/2mue) f Ieryrdr [ f(r, the~*P*dt. 
“sé ” 0 


The inversion of (2.11) can certainly be justified 
for the type of stresses that are likely to occur, 
since the stresses will be finite everywhere and 
will vanish beyond certain values of 7 and of ¢. 
With ¢(p, k) thus defined by (2.16), the question 
arises as to the types of media for which (2.12) 
converges and the subsequent inversion per- 
missible. We shall not enter here upon a dis- 
cussion of this question, and shall assume the 


(2.15) 


—oo 


(2.16) 





convergence of the integrals involved.® Equations 
(2.15) and (2.16) determine Z,(0, k, p)/Z,(0, k, p) 
as a function of the parameters k and p. It is 
from this quantity, taken in conjunction with 
the defining Eq. (2.8), that we wish to determine 
p(z) and u(z). 


5 The question of convergence of the integrals appearing 
in this and the following sections is bound up with the 
problem of free vibrations of the medium and with the 
condition that (2.10), (3.10) and (3.11) shall represent di- 
verging waves. 
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Following Langer, let us investigate the 
asymptotic behavior of Z, for large values of k. 
With Z=e/ "42, Eq. (2.8) becomes 


U+U?+(u/u)U+ p2(p/u)—k?=0. (2.17) 
Assume 
U=Z/Z=kU,+ Uy+(1/k) Uy 
+(1/k)?U_2---,® (2.18) 


substitute in (2.17) ana equate the coefficients of 
each power of k to zero. The coefficient of k is 
(U,;?—1), the root of interest to us being U; = —1. 
With this value for U; the other equations are: 


u/p= —2Uy, (2.19) 
Uy—2U_,— Ue + P?(p/n) =0, (2.20) 
U_,—2U_.=0, (2.21) 
U_2—2U_3+U_=0, ete. (2.22) 


If we now develop [k+Z,(0, k, p)/Z,(0, k, p)] 
into a descending power series in k, we can 
identify the various coefficients with the corre- 
sponding U’s in (2.18). The U’s are, in general, 
functions of the parameter p*. We can see, 
however, from (2.19) that Up» is independent of 
p*, since it equals a function of z only. It follows 
from (2.20) that 

U_,= U_wtp?U_n, (2.23) 
and that 
Uy—2U_w— Ue=0, (2.24) 2U_n=p/p. (2.25) 


From (2.21) we now infer that 


U g™= U_»+p?U_o, (2.26) 
Uu«~2U_a=0, (2.27) 
U_1—2U_2=0. (2.28) 

Similarly we infer from (2.22) that 
U_;= U_w»+ p? U_s,+ pt U_s, (2.29) 


* The asymptotic solutions of (2.8) for large values of p 
2 ‘tr . 
are exieS erm) has, These correspond to the solutions under 
the ray-concept, since, by (2.1), they represent waves 
travelling with a velocity (u/p)}. 
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Ua -3U + U_’?=0, (2.30) 
Uo ~3U «+230 «° U_,,=0, (2.31) 
2U_x»= U_3;’, etc. (2.32) 


Turning now to (2.19), we see that it determines 
(u/u) for z=0. When this equation’ is differ- 
entiated, it yields in conjunction with (2.24) 
d/dz(z/pu) at z=0. Similarly, one can obtain the 
higher derivatives of (i/u) from (2.27), (2.30), 
etc., each taken in conjunction with the deriva- 
tives of the preceding equations. In the same 
manner (2.28) determines d/dz(p/)o, (2.31) and 
(2.28) once differentiated determine d?/dz*(p/u)o, 
etc. In this way we can obtain a power series 
development of uw and p at the origin. 

It should be noticed that yo has to be deter- 
mined by an examination of the surface material, 
while (p/u)o is given by (2.25); also that we 
make use only of the equations starting with 
either U_;., or | oe Of the other equations, 
those which begin with U_in(n>1 ) are irrelevant, 
while an equation of the type (2.32) imposes a 
relationship between observable quantities. This 
feature is, however, a direct consequence of the 
differential equation in much the same way as 
the fact that the asymptotic development of 
Z,(0, k, p)/Z,(0, k, p) for large values of k starts 
with the term (—). 

If the applied stress, instead of being an 
arbitrary function of the time, were of a definite 
frequency p, the resulting motion would be 
periodic and of the same frequency. The coeffi- 
cients U;, in (2.18) would no longer be functions 
of the parameter ~* and the only result of the 
analysis would be a power series expansion for 
u(z). The arbitrariness of the variation in time 
of the applied stress, which amounts to a 
simultaneous application of pulses of all fre- 
quencies, is then the circumstance which makes 
possible the determination of p(z). 

To sum up, we have shown that when a 
radially symmetrical twisting stress is applied 
at the surface of a heterogeneous medium, the 
functions p(z) and y(z) can be determined 
uniquely from the observed surface motion by a 
process which involves only the operations of 
integration and the expansion in a power series. 
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3. CASE OF A DILATATIONAL DISTURBANCE 


Let us start out with a discussion of the formal solution of the inverse problem for the case 
when a radially symmetrical normal stress (p,.) is applied at the surface. The dilatation now being 
different from zero, the motion will depend on p, uw and X. One would expect therefore that an 
analysis of the resulting motion would yield all these three functions. 

In a cylindrical system of coordinates the radial, azimuthal and vertical displacements are by (1.8) 

u=dF/dr, v=(1/r)dF/00, w=dF/dz+w’. (3.1) 
If we limit ourselves again to solutions with cylindrical symmetry, we shall have 
v=0F/0=p.o=0. (3.2) 


Consider first periodic solutions of (1.10) and (1.11): 
F(z, r, t) =e'”'R(r) f(z), w'(z, r, t) =e!” 'O(r) ¢(z). (3.3) 


Substitute in (1.10), drop the time factor and divide by R(r)f(z), obtaining after a rearrangement 
of terms 


[peot(A+2u) (f/f) + 2uf/f]+(A+2u)(R/R+(1/r)R/R]+[(A+n)(¢/f) +u¢/f]O/R=0. (3.4) 

The dots stand for differentiation with respect to the independent variable in question. Denoting 
the terms in the first and last brackets of (3.4), respectively, by u(z) and v(z), we have 

u(z)/v(z)= —[R/R+(1/r)R/R](A+ 2p) /v(z) —O/R. (3.5) 


Upon differentiating now (3.5) with respect to z and dividing the result by d/dz[(A+2y)/v(z) ], 
one obtains an equation one side of which is a function of z only and the other side a function of r 
only. It follows that 


R/R+(1/r)R/R=—k, 2. R=Jo(kr). (3.6) 


Substituting this result into (3.4) we find that Q/R must be a constant, which we take to be unity, 
since we can absorb any constant factor in g(z) of (3.3). Thus 


Q(r) = Jo(kr). (3.7) 
With (3.6) and (3.7), Eqs. (1.10) and (1.11) become, after dropping the common factors, 


P of +(A+2u)(fF—kf)+(A+nu) o+u(2f+ ¢)=0, (3.8) 


b?o(f+¢)+(A+2u)(f—R+ 9) — Rue tM f— kf + ¢) +2u(f+ ¢) =0. (3.9) 


Upon eliminating one of the variables, one obtains a differential equation of the fourth order in 
the other variable. We shall, however, later find that when seeking the asymptotic solutions for 
large values of k the laborious process of elimination can be avoided; also that there exist such 
solutions which for fixed values of » and 2(#0) are like e~**. Let these elementary solutions be 
denoted by f(z) and ¢(z) without subscripts. We now seek solutions for F and w’ of the type 


F(z, r, t)= { eimap [ N(p, k)Jo(kr) f(z, k, p)dk, 
—_ 9 


w' (2, r= f cmap | N(p, k)Jo(kr) o(s, k, p)dk. 
— 0 
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Hence u=0dF/ar= — f eirtdp f N(p, k)kIi(kr)f(z, k, p)dk, (3.10) 
ait 0 
w=dF /az+w'= f eirtd p [ N(p, k)Jo(kr) f(z, k, p) + (2, k, p) dk, (3.11) 
-_ «9 
6=—°F+0w /az= f eirtd p f N(p, k)Jo(kr) (f¥— kf + oldk, (3.12) 
—c 0 


Piz =)b+2ndw/dz= f etd p f N(p, k)Jo(kr)[(A+2u) (f+ ¢) — AR dk, (3.13) 
oo 0 


per= ul (w/ar) + (Ou/d2)]= —u J eid p | N(p, k)kIi(kr) (f+ e)dk. (3.14) 
e co e/0 


The vanishing of ~., everywhere on the surface requires’ that 


(2f+ ¢)o=0. (3.15) 


The function N(p, k) has to be chosen so that (3.13), taken for z=0, shall equal P(r, t) — the applied 
stress. If, as is likely to happen in an actual case, P(r, t) is bounded and vanishes beyond certain 
values of r and of t, we can invert (3.13) obtaining 


[(A+2u) (f+ ¢) —ARF oN (p, k)/k= (1/277) f Jo(kr)rdr f P(r, t)e-‘P'dt, (3.16) 
0 —oo 


which defines N(p, k). Assuming now the convergence of (3.10) and (3.11) and that they can be 
inverted for the value z=0, we shall have 


f(0, k, p)N(p, k)= - (1/2n) { Ierrdr [ u(0, r, t)e~**dt, (3.17) 
0 —o 
C/O, k, p) + ¢(0, k, p) IN(p, k)/k= (1/2n) f Soler rar J w(0, 7, t)e*?*dt. (3.18) 
0 —co 
The last three equations determine uniquely 
[f+e)/fl and [(f+¢)/fhb (3.19) 


as functions of the parameters p and k. 

We are now confronted with a task similar to that of section 2, namely, to determine p(z), u(z) 
and X(z) from the two functions in (3.19). For that purpose we must find the asymptotic solutions 
of (3.8) and (3.9) for large values of k. As a first step let us prove that there exist solutions of the 
form e~**, Write Y=¢+f, W being the elementary z-factor for w. Substitute in (3.8) and (3.9), divide 
by uw and call 


P?p/u=a, w/u=B, d/u=7, (3.20) 
obtaining af +f—R(y+2)f+(1+yv+BW+f)=0, (3.21) 
a + (y+ 2)~ — ky +1)f — kW + (y +78) (W— Rf) +28) =0. (3.22) 


7 See Eq. (3.47). 
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Further let fsel* Uae YaeS* Sd d/dz(y/f)=(S—U)(W/f). (3.23) 


Substitute in (3.21) and (3.22), divide by f and differentiate once the first equation thus obtained. 
From the resulting three equations (¥/f) can be eliminated. One thus obtains two equations in U 
and S and their derivatives up to the second. If we now assume for U the development (2.18) and 
a similar development for S, we shall find that the coefficients of (k*) in the above equations, when 
set to equal zero, yield the relations: 


[U2 —(y+2) 0(y+2)S2—1] = — UiSi(y+1)?, (3.24) [U2 —(y+2)][(v+1)(S2— UiS:)]=9, (3.25) 


of which the non-trivial solutions are U;=S,; and U;= +1. We are interested, of course, in the pair 
U,= —1, S:= —1. It follows that, with Y= —kg, we are justified in writing 


f(z) =e**U(z), g(z) =e**S(z), (3.26) 
U(z) = Uo(z)+(1/k) U_u(z)+-:-, S(z) = So(z) + (1/k)Si(z)++>>. (3.27) 
The form (3.26) is essentially the same as (3.23), but the functions U and S are not identical with 


the previously defined functions. By (3.26), Eqs. (3.21) and the difference between (3.21) and 
(3.22) become 


aU+U+pU+k[ —2U—(y+1)S—B(U+S)]+k(y+1)(S—U)=0, (3.28) 


a(U—S) + U—(y+2)S+BU—(y¥+78+28)S+k[(y+3)(S— U)+(y+7B+8)(S—U)]=0. (3.29) 


On putting the coefficients of the powers of k equal to zero, there result the relations 


Ur= So, 

(y+1)(S_1— U4) —(y+3) H—28U=0, (3.30) 
(y+3)(S_1— Ua) +(y+7B+8)(S1— Ua) — (y+1) Uo— (¥+8 +8) Uo=0, (3.31) 
(y+1)(S_¢— U_s) —2U_,—(y+1)S_1—6(U_1+S_1) +aUo+ Hh+6Ur=0, (3.32)o, ; 


(y+3)(S_2— U_s) + (y+¥B +8) (S_2— U_2) +a(U4—S1)+ U—(y+2)S+6U4 
—(y+78+28)S1=0, (3.33)o,1 
(y+1)(S_3— U_3) —2U_2—(y+1)S_2—8(U_2+S_2) +aU_:+ U_,+8U_1=0, (3.34)o, 1, 2 
(y+3)(S_3— U_s) + (y+78+8)(S_3— U_s) +a(U_:—S_s) + U_2—(y+2)S_2+6U_s 
—(y+7B+28)S.=0, etc. (3.35)o,1.2 
The values of these functions at z=0 are restricted by (3.15), which is equivalent to 


k(U+S)=U. (3.36) 


Hence we have for z= 
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Uy=S=0 (3.37); U_,+S_1= U) (3.38); 


U_.+S_.= U_, (3.39)o, 1; U_3+S_3= Rs etc. (3.40)o, 1,2 


The U_; and S_; are in general functions of the parameter (*), since a occurs in the defining differ- 
ential Eqs. (3.30) to (3.35). We can prove, however, that 


UL; = U_iot+ (p’) U_at _—— + (p*)'U_ii (3.41) 


and similarly for the S_;; also that the values of U_; and S_; at z=0 are polynomials of degree 
(t—1) in (p*). 

We notice first that (S_;— U_;) and (S_,— U_,) in (3.31) can be eliminated with the aid of (3.30). 
There results an homogeneous second-order differential equation for Uy with coefficients which are 
independent of a. Hence Us is independent of p”. Turning now to (3.33), we see that we can substitute 
for (S_2— U_2) and (S_2— U_s) from (3.32). Also (U_,;—S_1), S_1, and S_; can be expressed in terms 
of U_; and U, by (3.30). One thus obtains an inhomogeneous second order differential equation 
for U_, with coefficients which are independent of p? and with a right-hand side which is linear in p’. 
Since the general solution of an inhomogeneous differential equation is /inear in the term on the 
right-hand side, it follows that U_,;= U_w+p?U_x. By (3.30) this result holds also for S_;. Similarly 
one can obtain from (3.35) and the preceding equations an inhomogeneous second order differential 
equation for U_, with coefficients which are independent of p*, but with a right-hand side which is 
quadratic in ~*. In general, the differential equation for U_, or S_, will be inhomogeneous and of 
the second order ; p? will occur only on the right-hand side and to the power of . Hence (3.41) is true. 

The initial values of U_;; and S_;; will be zero.* This can be established for U_,, and S_,, from the 
Eqs. (3.38) and (3.30), the latter taken for z=0. One can solve for U_; and S_; in terms of Up» and 
hence they cannot contain any (p*) term. Similar results for U_s2 and S_2 can be established from 
(3.39) and (3.32), and so on for the other functions. An equation like (3.32), taken for z =0, therefore 
stands for two equations, one independent of (p*) and the other multiplying (p*). This multiplicity 
is designated by the subscripts. 

It has been shown previously that it is possible to determine [(f+¢)/f]o and [(f+ ¢)/flo as 
functions of k and p. When taken in conjunction with (3.15) these determine uniquely (f/f) and 
(g/g), which in turn determine U/U and S/S. Remembering that Uy =.S)=0, we have 


U/U=[Uo+(1/k) Ut: ++ J/((1/k) Ut (1/k)?U_2+ +++ J=kLitLo+(1/k)Lit---, (3.42) 
S/S=[So+(1/k)Sat- ++ V/[(1/k)S1+(1/k)2S_o+ +++ J=kRM+My+(1/k)Myt---, (3.43) 
implying that 
(Uo—L,U_s) + (1/k)[U_1—(Li1U_2+LoU_1)] 
+(1/k)?{ U_e—(L,U_3+LoU_2+L_1~U_:)]+-::=0, (3.44) 


(So— M,S_:) + (1/k)[S_1— (M1 S_2+ MoS_1)] 


+ (1/k)*{S_2—(M,S_3+ MoS_2+M_,S-1)]+--+=0. (3.45) 


The L; and M, are determined experimentally as functions of p. They are defined in terms of the 
U’s and S's by putting every coefficient of (1/k)" in (3.44) and (3.45) equal to zero. Now since 
Uy and U_, are independent of (p?), L; is independent of (p?). The (1/k) term in (3.44) shows that 


‘It is evident from the system of Eqs. (3.30) to (3.35) that U_ii =S_i for all values of z. 
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Ly may be linear in (p*); generally the L_, is a polynomial of degree (k+1) in (p*). The same is 
evidently true for the M_,. The coefficients of (1/k)”" in (3.44) and (3.45), when put equal to zero, 
stand therefore for (n+1) equations corresponding to the required vanishing of the individual 
coefficients of the powers of (p”). 

It is now possible to determine a, 8 and y as functions of z from the system of equations consisting 
of (3.44) and (3.45), and the groups of equations starting with (3.37) and (3.30), respectively, the 
latter taken for z=0. The procedure is indicated in Table I where we have written on the left side 


TABLE I. List of unknowns and groups of equations from which they can be determined. 








Sn Sia Gee (3.30), (3.38), (k°) terms in (3.44) and (3.45) 

U_20, U_o1, S_20, S-21, U_to, (3.30), (3.31), (3.32)o, 1, (3.39)o, 1, (1/k) terms in (3.44)o, 1, (3.45)o, 1 
U_n, S_10, Uo, 7’: B 

U_s0, U_s, U_se, S-s0, Ss, S_s2, (3.30), (3.31), (3.32)o, 1, (3.33)o, 1, (3.34)o, 1, 2» (3.40)o, 1, 2 

U_20, U_a, U_22, S_20, S_21, (1/k)? terms in (3.44)o, 1,2 and (3.45)o, a3 


ses ee . 


U_10, U_11, S10, Uo, y, B, & 








groups of unknowns and on the right side the corresponding groups of equations from which they 
can be determined. Since the set of equations starting with (3.30) holds for all values of z, we can 
differentiate them and take the resulting equations for z=0. Such equations are signified by dots 
over the numbers of the equations. One thus finds, for example, that 


y= —1-—JL), B= (1/2L;)[Lo(Li+1)+Mo(Li—1)?], y= (4)[Lo(Li—1) — Mo(Li—1)*], etc. 


As in the case of non-dilatational disturbance one meets here with the circumstance that certain 
relations exist between observable quantities. Thus we must have L;—1=L;/M,;° also the number 
of the next group of equations which determines ¥, 8 and & exceeds the number of unknowns by 
one. This phenomenon is, again, an intrinsic characteristic of the differential equations and will 
lead to identities in an actual case. 

The above theory can easily be extended to the case where simultaneously with the normal 
stress (p..) P(r, t) there is applied a radially symmetrical radial stress (p.,) —Q(r, ¢). Instead of Eq. 
(3.15) we sha'l then have, by (3.14), 


— f eirtdp f N(p, k)kIi(kr)(2f+ g)odk=O(r, t), (3.46) 


which when inverted, yields the relation 


N(p, k)(2f+ y)o= — (1/2mp0) f Ti(kr)rdr | Q(r, te-‘”*dt. (3.47) 


Eq. (3.47) in conjunction with (3.16), (3.17) and (3.18) determines uniquely (f/f)o, (¥/W)o and (W/f)o, 
as functions of k and p. These in turn determine (U/U)o, (S/S)o and (S/U)p. Since Up and So need 


not be zero now we shall have, in general, 


U/U=Mot+(1/k)Mit:--, S/S=No+(1/k)Nat-::, S/U=14+(1/k)Art:::, (3.48) 


the M’s, N’s and A’s being defined by putting the coefficients of the powers of (1/) in 


* The relation y = —1—LZy, is of the same type, since y =(A/u)o can be determined independently. 
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(Uo— UsMo) + (1/k) [U1 — (Uo M_1+ U_1Mo)] 


+ (1/k)*{ U_2—(UoM_2+ U_:M_1+ U_2Mo)]+->-, 


(So—SoNo)+(1/k)[S—1— (SoN_1+S_1No)]+ (1/k)?[S_2— (SoN_-2+S_1N1+S_2No)]+-->, 


and in 


(1/k)[.S_1— (U_1+ UsA_1) J+ (1/k)*(S_2— (U_2+ U_,A_1+ UpA-_2)]+:>>, 
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(3.49) 


(3.50) 


(3.51) 


equal to zero. Eq. (3.41) will still be true, but U_;;(=S_;;) need not be zero at the surface. 

With the M’s, N’s and A’s known as functions of (p?), the determination of \, w-and p can be 
carried out as in the previous section, with the only difference that Eqs. (3.49) and (3.50) now take 
the place of (3.44) and (3.45), while Eq. (3.51) replaces the system of equations starting with (3.37). 


One thus finds 


B= (3) [(y+1)A-1— (y¥ +3) Mo], 


B= (y+1)A-2—B(A-10+8+ Mo) + (y+2)A-1w.Me— M_w— (y+2)N_.0, — = (y+1)A-a+(y¥+3)M_y 
= [2(y+2)(N_10— M_10-+ MoA—10) + (y+1)(8+8Mo) | 


+ (3)(v¥+3)(v¥+1)B(A-10— Mo) J/(Mo—A-10), ete. 


It follows that when the (radially symmetrical) normal and radial stresses applied at the surface 
are known, it is possible to determine from the observed surface motion the variation with depth 


of p, uw and X. 


4. CONCLUSIONS 


We can sum up our results as follows: The 
motion of a heterogeneous half-space resolves 
itself, as in a homogeneous medium, into two 
types: one—unaccompanied by dilatation, and 
one which is accompanied by dilatation. If the 
motion is to have cylindrical symmetry, the 
stresses applied at the surface must also have 
cylindrical symmetry and, furthermore, in the 
case of a non-dilatational disturbance the only 
possible stress is the azimuthal (f.9), while in 
the case of a dilatational disturbance, the possible 
stresses are the normal (p,,) and the radial (p.,). 
In the former case the motion is horizontal and 





in the azimuthal direction (v), in the second 
case the motion takes place in planes through 
the vertical, the azimuthal displacement being 
zero. If normal, radial and azimuthal stresses, 
of known space and time distributions, are 
applied simultaneously at the surface and the 
subsequent surface motion is observed, one can 
determine p(z) and y»(z) from v alone and p(z), 
u(z) and X(z) from uw and w alone. 

It is a pleasure to acknowledge my indebted- 
ness to Professor L. B. Slichter who suggested 
the problem and directed the research. I also 
wish to express my appreciation to Professor 
R. E. Langer for having criticised the paper in 
manuscript. 
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Study of Emergence Angle and Propagation Paths of Seismic Waves 
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Travel-time curves are given for elastic waves through 
the earth in a region in which the velocity increases 
continuously with depth. These time-distance curves have 
been approximated by an equation of the form X =a7*+bT 
from which the velocity depth relation has been deduced. 
A new method for measurement of the emergence angle 
has been used and the values obtained agree reasonably 
well with those deduced from theoretical treatment of the 
time-distance curve. Formulas for the travel-time between 


any two points in the medium under consideration are 
derived. These formulas agree closely with two sets of 
direct measurements. In one of these the seismograph was 
placed at various distances directly beneath the explosion. 
In the second case it was placed at a fixed depth beneath 
the surface and the explosion was located on the surface 
at various distances. The approximate depth to bed rock 
is obtained from the time-distance curves by use of the 
formulas mentioned above. 





HE theory of the propagation of elastic 

waves through the earth in a region in 
which the velocity increases with depth in such 
a fashion that the observed time-distance curve 
may be represented by the equation 


X =aT*+6T (1) 


has been treated by Ewing and Leet.' The 
present paper includes ari extension of the 
theory and an experimental investigation of 
many points in connection therewith. The data 
used were taken near Green Pond, Northampton 
County, Pennsylvania. The material at the 
surface was a fine grained alluvial deposit, which 
was underlain by limestone. The surface of the 
limestone was probably quite irregular. 


TIME-DISTANCE DATA AND VELOCITY-DEPTH 
DETERMINATION 


The time-distance data are shown in Fig. 1. 
The curved lines drawn represent Eq. (1) with 


a= 29,690 ft./sec.?, b=752.9 ft./sec. 


The values of these constants were obtained by 
a least square solution which included all points 
representing wave paths lying entirely in the 
upper layer. By use of the following equations 


7=2aT+b, (2) 
b/sin 09=v/sin 6= V, (3) 
P=6*/8ra(sinh g—q),' (4) 


g=2 arc cosh V/d 


‘Ewing and Leet, Trans. A.I.M.E. Geophysical Pros- 


pecting (1932). 


the velocity depth relation is determined. Here V 
is the velocity at the midpoint of the path of 
length X and maximum penetration P, and 4 is 
the emergence angle for this path. Fig. 2 illus- 
trates the meaning of each of these symbols. The 
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Fic. 1. Time-distance curves and approximate section at 
Green Pond. 
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Fic. 2. Diagram of wave path. 
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Fic. 3. Velocity-depth relation (Green Pond). 


velocity-depth curve thus obtained is shown in 
Fig. 3. 


EMERGENCE ANGLE 


A new method was used to measure the 
emergence angle 4) of the waves. The seismo- 
graph was of the type technically known as a 
geophone. It was originally designed to measure 
only vertical motions but after an alteration of 
the spring suspension was made it would operate 
when tilted at any angle between the vertical 
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Fic. 4. Seismograms illustrating method used to de- 
termine emergence angle. Angles given represent tilt of 
seismograph from vertical position. A, low sensitivity; B, 
high sensitivity; C, timing line. 











AND A. P. CRARY 
160} Observed © 
Computed — 
- (e0r 
+ 
v 
z 
© 80 
2 
=) 
40Fr 











4 iL i 
40° 30° 20° 10° 
Emergence Angles 


Fic. 5. Relation of emergence angle to distance traversed. 


and the horizontal. The tilting produced some 
change in the constants of the instrument, but 
this change was of no importance in the inves- 
tigation. The idea of the experiment was to tilt 
the geophone toward the shot until its axis 
became perpendicular to the direction of the 
first motion produced by the shot. Successive 
tests were made to determine the angle of tilt 
at which the direction of first motion on the 
seismogram showed a reversal. (This angle could 
easily be determined to one or two degrees in 
the range of distances studied.) Fig. 4 shows the 
set of observations taken at 5 feet. In Fig. 5 
emergence angle 4, is plotted against the distance 
traversed by the wave along the surface. The 
curve drawn represents the theoretical values of 
the angle as determined from Eq. (3) using the 
values of a and } obtained from the time-distance 
curves. The observed values, indicated by the 
points, agree with the calculated values within 
the limits of experimental error. It was not con- 
sidered necessary to correct for the difference 
between true and apparent emergence angles.’ 


TRAVEL-TIME BETWEEN ANy Two PoINTs ON A 
PATH 
From Eq. (4) the relation 
dy/dv = (v? —b*)!/27a, 


2 E. Wiechert, Gott. Nachr. 1907, No. I. 
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PROPAGATION OF 


where v is the velocity at the depth y, is obtained. 
Eq. (2) yields 


dy/dx = +(B csc? 0)—v")!/2, (6) 


which may be combined with (5) to give the 
relation: 


x= +(1, 2ra) f v[_(v? — b?) /(b? csc? 09 —v*) ]idv 
b 


= (b? cot? 0)/4ra)La—sin acosal, (7) 


where sin a= -ttan 69(v?/b®?—1). The positive 
sign should be taken if the wave has not reached 
its maximum depth. 

The travel-time /= {ds/v becomes, upon sub- 
stitution from (5) and (6) 


t= +(b/2za sin wo) f (1/v) 
vp 


x L(x? — b?) /(B? csc? 0) —v*) Jide 
= (b/27a_.sin 9) 
XLa—sin 4) tan“'(tan a/sin 4) ]. (8) 


For the limiting case, 6>=0, the path is ver- 
tical and (7) and (8) become 


x=0 (7’) 
t=(b/27a)[(v?/b?—1)!—tan-"(v?/b?—-1)!].  (8’) 


Formulas (7) and (8) make it possible to cal- 
culate the time required to traverse any part of 
any path. Two sets of data were taken to which 
these formulae may be applied. In the first case 
records were taken with the geophone buried at 
a series of depths up to about 20 ft., the shot in 
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Fic. 6. Travel time curve for waves travelling vertically. 
Fic. 7. Travel time curve for waves travelling from surface 
to depth of 18’ 6”. 
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each case being at the surface directly above. In 
the second case the geophone remained at a 
depth of 18.5 ft. while shots were fired on the 
surface at various distances from the mouth of 
the hole. Figs. 6 and 7 compare the travel-time 
obtained from these observations with those 
computed from the above equations using the 
values of a and } as computed from the time- 
distance curves. 


WAVES REFLECTED OR REFRACTED AT 
DISCONTINUITIES 


By use of the formulas developed in the pre- 
ceding section it is possible to calculate time- 
distance curves for waves reflected or refracted 
at the surface of a layer of rocks covered by 
matetial in which the velocity increases with 
depth in the manner discussed above. For a 
given depth the velocity v is given by Eq. (4). 
From Egs. (7) and (8) it is possible to calculate 
time-distance curves for waves reflected at this 
depth by assuming a series of values for 0o. If the 
reflecting surface is horizontal the time and the 
distance will be just double that calculated from 
these equations. If the 1eflecting surface is sloping 
it is necessary to make separate calculations for 
the incident and the reflected waves. 

In the case of waves which are refracted along 
the surface of a buried rock the angle of emer- 
gence is known from the velocity of the refracting 
layer. If, then, a depth is assumed the horizontal 
distance and the time of the wave in the upper 
layer can be computed from the value of the 
emergence angle and the velocity in the upper 
layer corresponding to the depth assumed, as in 
the case of reflected waves. The remainder of the 
distance traversed is in the higher speed material 
and the total time is the sum of the time spent 
there and the time of the incident and emergent 
rays in the upper layer.* For waves refracted 
along a sloping layer, the computations become 
increasingly more difficult as the emergence angle 
depends not only on the velocity of the high 
speed layer but also upon its depth. 

Although the primary purpose of the present 
work was to study waves propagated in the 
upper layer, the data represented in Fig. 1 afford 
a good illustration of refracted waves. The 


3 Rutherford, Am. Geophys. Union. Trans., p. 292, 1933. 
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straight lines drawn in Profiles 3 and 4 represent 
an assumed horizontal layer which is 100 ft. 
beneath the surface and has a velocity of 17,000 


ft./sec. Because of irregularities in the surface of 
the rock the observed points deviate considerably 
from the lines, but the depth assumed is certainly 
accurate to 5 or 10 percent. Profiles 1 and 2 
indicate a limestone surface sloping up sharply 
to the east. By the aid of Fig. 6 the depth at 
the east end may be estimated at approximately 
20 ft. In the region between the 500 and 600 ft. 
stations there are systematic deviations marked 
by dotted lines on the time-distance curve. These 
indicate a fairly sharp depression in the surface 
of the limestone between these stations. The 
deviations of the points in general are an indica- 
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CRARY 


tion of the roughness of the limestone surface. 
An approximate cross section of this region is 
shown in Fig. 1. 

Professor B. L. Miller of the Department of 
Geology, of Lehigh University, who recom- 
mended the location which was chosen for these 
experiments, has suggested that the structure of 
the limestone revealed by these results represents 
the buried valley of a pre-glacial stream. 

We are indebted to the Geophysical Research 
Corporation for the seismic instruments, to the 
Trojan Powder Company for the explosives, to 
Mr. H. S. Snyder upon whose property the ex- 
periments were carried out, and to Mr. A. M. 
Thorne, Jr., for assistance with the experimental 
observations. 
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